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The  NATO  Research  Technology  Organization  (RTO)  Advanced  Vehicle  Technology 
(AVT)  Panel  established  Working  Group  10  (WG  10)  “Technologies  for  Propelled  Hy¬ 
personic  flight”  in  1998.  The  Terms  of  Reference  included  three  Subgroups:  1)  aerospike 
plug  nozzle,  2)  RAM/SCRAM  jet  conflgurations,  and  3)  Computational  Fluid  Dynamics 
(CFD)  validation.  An  overview  of  Subgroup  3  (SG  3)  is  presented  in  this  paper.  The 
SG  3  participants  deflned  six  topical  areas  for  which  validation  of  CFD  methodologies  was 
deemed  essential  for  eflfective  analysis  and  design  of  propelled  hypersonic  vehicles.  The 
topics  are  boundary  layer  instability  and  transition,  real  gas  flows,  laminar  hypersonic 
viscous-inviscid  interactions,  shock-shock  interactions,  shock  wave-turbulent  boundary 
layer  interaction,  and  base  flows  with  and  without  plume  interaction.  A  set  of  experi¬ 
ments  were  selected  by  Subgroup  3  participants  for  each  topic  and  computations  were 
solicited.  The  experiments  are  described  in  this  paper.  The  results  of  the  computations 
are  presented  in  separate  papers. 


Introduction 

The  objective  of  this  paper  is  to  document  a  set  of 
experiments  for  validation  of  CFD  capability  for  spe¬ 
cific  flow  phenomena  relevant  to  hypersonic  flight.  The 
topic  areas  were  established  by  NATO  RTO  Work¬ 
ing  Group  10  Subgroup  3  at  its  first  meeting  in  De¬ 
cember  1998  in  Chalais-Meudon,  France,  and  final¬ 
ized  at  its  meetings  in  June  1999  in  Norfolk,  Vir¬ 
ginia,  and  November  1999  in  Lampoldshausen,  Ger¬ 
many.  The  topic  areas  and  team  leaders  are  listed 
in  Table  1.  In  each  area,  specific  experiments  were 
selected  by  the  Subgroup  for  CFD  validation.  The 
team  leaders  for  each  topic  prepared  the  sections  in 
this  paper  and  solicited  participants  to  contribute  one 
or  more  computations.  A  website  is  available  at: 
http :  /  /  WWW.  engr  .rutgers .  edu/  ~  wgl  0 . 

The  following  sections  contain  a  description  of 
the  selected  experiments.  In  some  cases,  nominal 
freestream  test  conditions  are  provided  and  do  not  nec¬ 
essarily  represent  the  exact  freestream  test  conditions 
for  which  experimental  data  was  obtained.  The  indi¬ 
vidual  researchers  should  be  contacted  in  all  cases  for 
the  exact  freestream  test  conditions. 

A  summary  of  the  computations  for  each  of  the  topic 
areas  are  presented  in  a  series  of  AIAA  papers  at  the 
AIAA  39th  and  40th  Aerospace  Sciences  Meetings. 
The  list  of  papers  is  shown  in  Table  2. 

A  selected  set  of  computations  are  described  in  de¬ 
tail  in  a  series  of  papers  at  the  AIAA  40th  Aerospace 
Sciences  Meeting.  The  list  of  papers  is  shown  in  Ta¬ 
ble  3. 
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1.  BOUNDARY  LAYER  INSTABILITY 
AND  TRANSITION 
TEAM  LEADER:  S.  SCHNEIDER 

Prepared  by 

S.  Schneider 

General  Introduction 

Laminar-turbulent  transition  has  a  significant  in¬ 
fluence  on  hypersonic  heat  transfer,  boundary-layer 
separation,  and  drag.  Transition  is  affected  by  many 
factors,  including  freestream  noise,  roughness,  wavi¬ 
ness,  geometry,  Mach  number,  and  wall  temperature. 
Possible  mechanisms  include  the  first  and  second-mode 
instabilities,  the  3D  crossflow  instability,  the  concave- 
wall  Gdrtler  instability,  and  streamwise  vorticity.  Saric 
et  aO  present  a  recent  review.  Experimental  data  that 
measures  only  transition  location  is  difficult  to  analyze 
since  the  mechanism  is  often  completely  unknown;  the 
systematic  uncertainty  cannot  in  general  be  resolved. 
Gorrelations  of  transition  which  ignore  the  transition 
mechanisms  can  only  be  valid  for  the  conditions  of  the 
pre-existing  correlated  data.  Attempts  to  correlate  a 
wide  range  of  transition  data  using  algebraic  formulas 
suffer  from  scatter  that  spans  an  order  of  magnitude 
(see  Schneider^). 

Reliable  predictive  methods  will  have  to  be  based  on 
simulations  of  the  transition  mechanisms.  The  best 
developed  of  these  mechanism-based  methods  is  the 
method,  which  compares  the  integrated  growth 
of  the  dominant  instability  waves  to  empirical  transi¬ 
tion  data.  More  advanced  mechanism-based  methods 
include  Direct  Numerical  Simulations  or  use  of  the 
Parabolized  Stability  Equations.  A  necessary  step  in 
the  further  development  of  these  methods  is  accurate 
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Topic  Title 
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2  Real  Gas  Flows 

3  Laminar  Viscous-Inviscid  Interactions 

4  Shock-Shock  Interactions 

5  Shock  Wave- Turbulent  Boundary  Layer  Interactions 

6  Base  Flows  with  and  without  Plume  Interactions 


Team  Leader (s) 

Steven  Schneider  (Purdue  Universty) 

Graham  Candler  (University  of  Minnesota) 
Michael  Holden  (CUBRC) 

Steven  Walker  and  John  Schmisseur  (AFOSR) 
Doyle  Knight  (Rutgers  University) 

Peter  Bakker  (Technical  University  of  Delft) 
Phillipe  Reijasse  (ONERA) 


Table  2  Topics  and  Papers 


Topic  Author(s) 

ov  Jean  Muylaert 

1  Steven  Schneider 

2  Graham  Candler 

3  Michael  Holden  and 
Timothy  Wadhams 


3  John  Harvey,  Michael  Holden 
and  Timothy  Wadhams 


3  Michael  Holden  and  John  Harvey 


4  Steven  Walker  and 
John  Schmisseur 

5  Doyle  Knight,  Hong  Yan, 

Argyris  Panaras  and 
Alexander  Zheltovodov 

6  Peter  Bakker,  Willem  Bannink 
and  Phillipe  Reijasse 

ov  Overview  of  RTO  Working  Group  10 
*  There  are  two  papers 


Title 

RTO  Working  Group  10:  Technologies  for 
Propelled  Hypersonic  Flight 
Hypersonic  Laminar  Instability  on  Round 
Cones  Near  Zero  Angle  of  Attack 
CFD  Validation  for  Hypersonic  Flight: 

Real  Gas  Flows 

Code  Validation  Study  of  Laminar 
Shock/Boundary  Layer  and  Shock/Shock 
Interactions  in  Hypersonic  Flow 
Part  A:  Experimental  Measurements 
Code  Validation  Study  of  Laminar 
Shock/Boundary  Layer  and  Shock/Shock 
Interactions  in  Hypersonic  Elow 
Part  B:  Comparison  with  Navier-Stokes 
and  DSMC  Solutions 

Comparisons  Between  DSMC  and  Navier-Stokes 
Solutions  on  Measurements  in  Regions  of 
Laminar  Shock  Wave  Boundary  Layer 
Interaction  in  Hypersonic  Elow 
CED  Validation  of  Shock  Shock  Interaction 
Elowfields 

RTO  WG  10:  CED  Validation  for  Shock 
Wave  Turbulent  Boundary  Layer  Interactions 

Base  Elows  With  and  Without  Plume 
Interactions 


Paper  No. 
AIAA  2002-0432 

AIAA  2001-0206 

AIAA  2002-0434 

AIAA  2001-1031* 


AIAA  2001-1031* 


AIAA  2002-0435 

AIAA  2002-0436 
AIAA  2002-0437 

AIAA  2002-0438 


Topic 

2 

4 
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Table  3  Topics  and  Papers 


Author  (s) 

Graham  Candler,  1.  Nompelis, 
M.  C.  Druguet,  1.  Boyd  and 
W.-L.  Wang 
Dominic  D’Ambrosio 

Erederic  Thivet 


Title 

CED  Validation  for  Hypersonic  Elight 
Hypersonic  Double-Cone  Elow  Simulations 

Numerical  Prediction  of  Shock  /  Shock 
Interactions  in  Hypersonic  Elow 
Lessons  Learned  Erom  RANS  Simulations 
of  3D  Shock  Wave  Boundary  Layer  Interactions 


Paper  No. 
AIAA  2002-0581 

AIAA  2002-0582 
AIAA  2002-0583 
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computation  of  the  actual  instability-wave  growth, 
which  is  not  trivial.  For  hypersonic  flight  vehicles  with 
reasonably  smooth  walls,  transition  is  likely  to  occur 
via  the  growth  of  linear-instability  waves,  so  methods 
such  as  are  an  appropriate  approach  to  prediction. 

Most  ground-test  data  are  contaminated  by  the  high 
levels  of  noise  which  normally  radiate  from  the  turbu¬ 
lent  boundary  layers  on  the  nozzle  walls.®  Even  in 
these  noisy  conventional  tunnels,  there  are  few  mea¬ 
surements  of  the  instability  mechanisms  causing  tran¬ 
sition.  Of  these  instability  measurements,  few  have 
been  carried  out  with  calibrated  instruments  and  well- 
documented  conditions,  so  that  a  fairly  reliable  com¬ 
parison  to  computations  can  be  carried  out. 

The  team  feels  that  the  following  elements  should  be 
present  in  code-validation  experiments  for  hypersonic 
transition: 

•  Detailed  and  reliable  measurements  of  the  transi¬ 
tion  mechanisms 

•  Accurate  knowledge  of  the  mean  flow 

This  includes  repeated  checks  of  boundary  layer 
symmetry,  instrumentation  calibrations,  tunnel 
flow  nonuniformity  effects,  repeatability,  etc. 
Moreover,  close  agreement  between  computation 
and  experiment  for  mean  flow  characteristics  is  es¬ 
sential.  This  requires  repeated  tunnel  entries  and 
cooperation  between  experiments  and  computa¬ 
tions. 

•  Accurate  measurements  of  the  fluctuating  flow- 
field 

Comparison  with  linear  instability  theory  requires 
calibrated  measurements  of  the  fluctuations  at  a 
known  position  in  the  eigenfunction.  It  seems 
possible  that  these  might  be  carried  out  in  con¬ 
ventional  tunnels  using  ensemble  averaging  of  con¬ 
trolled  disturbances. 

•  Experiments  in  quiet  tunnels 

Nonlinear  secondary  breakdown  effects  are  depen¬ 
dent  on  small  fluctuations  combining  with  the 
primary  instability.  It  seems  doubtful  that  these 
can  be  repeatably  and  successfully  studied  except 
in  quiet  tunnels,  since  even  controlled  secondary 
disturbances  may  be  swamped  by  tunnel  noise. 
Eor  similar  reasons,  receptivity  experiments  may 
also  require  quiet  flow. 

These  specifications  are  updated  from  Reshotko.^ 
This  earlier  reference  gives  additional  suggestions,  and 
shows  the  long-term  importance  of  a  coordinated  ap¬ 
proach.  Unfortunately,  these  specifications  are  ex¬ 
tremely  difficult  to  meet  at  hypersonic  speeds.  Al¬ 
though  new  experimental  efforts  are  underway  at  Pur¬ 
due  and  at  the  Institute  of  Theoretical  and  Applied 


Mechanics  (ITAM),  Novosibirsk,  Russia,  all  existing 
data  fall  short. 

Three  experimental  datasets  are  described  below. 
A  detailed  analysis  of  the  datasets  and  prior  com¬ 
putations  is  presented  in  Schneider.®  While  all  three 
datasets  contain  significant  flaws  and  limitations,  they 
are  the  best  presently  available,  in  the  opinion  of  the 
team.  Accurate  computation  of  the  wave  growth  and 
eigenfunctions  for  these  cases  would  be  a  significant 
step  towards  the  further  development  of  mechanism- 
based  methods.  Comparisons  between  the  experi¬ 
ments  and  the  results  of  multiple  independent  com¬ 
putations  would  shed  light  on  the  strong  and  weak 
elements  in  all,  and  aid  in  the  development  of  new 
work. 

Because  of  the  limited  resources  presently  available, 
attention  is  focused  on  axisymmetric  geometries.  The 
computational  resources  required  are  much  less  than 
in  the  3D  case.  All  of  the  datasets  are  for  cold  hyper¬ 
sonic  flow  without  chemistry,  because  detailed  mea¬ 
surements  of  the  mechanisms  of  hot  hypersonic  transi¬ 
tion  do  not  exist.  Although  the  data  shown  below  are 
primarily  measurements  of  the  cold-flow  second-mode 
instability,  other  mechanisms  do  of  course  remain  im¬ 
portant  (e.g.,  Hornung®). 

Dataset  No.  1:  Sharp  Cone  at  Mach  8 

Prepared  by 

S.  Schneider 

Introduction 

Stetson  et  ®  carried  out  calibrated  measurements 
of  instability  wave  growth  on  several  conical  models  in 
AEDC  Tunnel  B  at  Mach  8,  during  the  late  1970’s  and 
early  1980’s.  Detailed  hot-wire  measurements  were 
carried  out  in  this  expensive  production  tunnel  by 
J.  Donaldson,  under  the  direction  of  K.  Stetson.  The 
experiments  focused  on  the  hypersonic  second-mode 
instability,  which  is  likely  to  be  dominant  on  smooth 
convex  nearly-symmetric  geometries  with  small  cross- 
flow.  Second-harmonic  amplification  was  observed 
downstream,  a  good  indication  of  nonlinear  effects.  A 
comparison  of  the  second-harmonic  amplification  rate 
to  a  nonlinear  computation  would  be  interesting.  A 
number  of  workers  have  already  compared  these  results 
to  computations;  however,  none  have  taken  advantage 
of  the  substantial  amount  of  tabulated  data  that  is  not 
presented  in  the  conference  papers.  Although  these  are 
conventional-tunnel  measurements  with  high  ambient 
noise,  the  good-quality  detailed  measurements  appear 
to  deserve  further  computational  comparisons.® 

Brief  Description  of  Experiments 

The  model  was  a  7  degree  half-angle  cone  with  a 
40  inch  (1.016  m)  length  and  a  sharp  0.0015  inch 
(38  micron)  nose  radius.  The  cone  angle  of  attack 
was  zero  (within  the  accuracy  with  which  this  could 
be  set  and  measured).  The  model  was  in  thermal  equi- 
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librium  but  heat  conduction  within  the  model  may  not 
be  negligible. 

Most  of  the  measurements  were  carried  out  at  a 
freestream  Mach  number  of  7.95,  with  a  total  pres¬ 
sure  of  225  psia  (1.55  MPa)  and  a  total  temperature 
of  1310  R  (728  K).  Additional  data  was  recorded 
at  a  lower  total  pressure  (to  determine  the  neutral- 
instability  region)  and  at  a  higher  pressure  (to  record 
data  into  transition  onset).  Thus,  transition-onset 
data  should  be  available. 

Experimental  data  included  surface  pressure  taps 
(24)  and  thermocouples  (32),  mean-flow  profiles  us¬ 
ing  pitot-pressure  and  total  temperature  probes  and 
instability-wave  spectra  using  calibrated  constant- 
current  hot-wire  anemometry.  Hot-wire  measurements 
were  recorded  at  one-inch  (2.54-cm)  intervals  between 
10  and  37  inches  (25.4  and  94.0  cm),  using  a  single 
hot-wire  with  a  single  calibration. 

The  experimental  data  is  available  on  the  RTO 
WG  10  SG  3  website  (see  above).  See  Schneider®  for 
a  detailed  discussion. 

Dataset  No.  2:  Blunt  Cone  at  Mach  8 

Prepared  by 

S.  Schneider 

Introduction 

In  addition  to  the  sharp  cone  measurements  in 
Dataset  No.  1,  Stetson  et  al^’^  carried  out  additional 
calibrated  measurements  of  instability  wave  growth 
on  several  blunt  cone  models  in  AEDC  Tunnel  B  at 
Mach  8.  The  apparatus  was  almost  identical.  This 
conventional  wind  tunnel  is  well  characterized. The 
blunt-cone  experiments  examined  the  effects  of  the  en¬ 
tropy  layer  on  the  hypersonic  second-mode  instability. 
Herbert  et  al,^^  Kufner  et  and  others  have  al¬ 
ready  compared  some  of  these  results  to  computations. 
Eig.  2  in  Kufner  et  al  summarizes  the  comparisons. 
All  of  the  computations  show  a  growth  rate  about  1.5 
times  the  experimental  value,  at  a  frequency  roughly 
10%  below  the  experimental  value.  However,  these 
comparisons  were  carried  out  at  a  station  175  nose 
radii  downstream  of  the  tip.  Eig.  7b  in  Stetson  et 
al^  paper  shows  substantial  amplification  of  second- 
harmonic  frequencies  at  this  station.  Stetson  believes 
that  nonlinearity  is  substantial  at  this  station,  there¬ 
fore  a  growth-rate  comparison  to  linear  theory  has  du¬ 
bious  validity.  Many  workers  have  selected  this  station 
for  comparison  to  linear  theory,  apparently  through  in¬ 
sufficient  communication  with  Stetson  et  al. 

The  mean  flow  was  measured  for  this  experiment, 
but  existing  comparisons  have  uncovered  some  diffi¬ 
culties.  Eig.  5  in  Herbert  et  al^^  shows  that  an 
adiabatic- wall  computation  overpredicts  the  wall  tem¬ 
perature  by  roughly  20%.  This  is  thought  to  be  due  to 
heat  transfer  to  the  water-cooled  sting  support.  While 
the  measurements  were  carried  out  after  thermal  equi¬ 


librium  was  achieved,  heat  transfer  within  the  model 
apparently  cannot  be  neglected. 

Boundary-layer  profiles  were  also  measured  using 
pitot  and  total-temperature  probes.  Unfortunately, 
repeated  runs  reveal  discrepancies  of  10-15%  under 
nominally  the  same  conditions  {e.g.,  Eig.  6  in  Her- 
bert^^).  Also,  the  data  probably  suffer  from  probe 
interference  effects.® 

The  data  show  growth  of  the  first  and  second-mode 
instability  waves  on  a  symmetric  blunt  cone  under 
cold-flow  conditions.  The  second-mode  instability  is 
dominant.  Although  amplification  was  only  observed 
downstream  of  the  nominal  swallowing  length,  strong 
bluntness  effects  were  apparent.  Some  entropy-layer 
fluctuations  were  also  apparent,  although  these  could 
not  be  resolved  clearly  and  no  quantitative  compar¬ 
isons  have  been  attempted.  Second-harmonic  am¬ 
plification  of  the  second-mode  waves  was  observed 
downstream,  a  good  indication  of  nonlinear  effects.  A 
comparison  of  the  second-mode  amplification  rate  to  a 
nonlinear  computation  would  be  interesting. 

Description  of  Experiment 

The  model  was  a  7-degree  half-angle  cone  with  a  40- 
inch  (1.016-m)  length  and  a  0.15-inch  (3.81  mm)  nose 
radius.  The  cone  angle  of  attack  was  zero,  to  within 
the  accuracy  with  which  this  could  be  set  and  mea¬ 
sured.  The  model  was  in  thermal  equilibrium  but  heat 
conduction  within  the  model  may  not  be  negligible. 

Most  of  the  measurements  were  carried  out  at  a 
freestream  Mach  number  of  7.99,  with  a  total  pres¬ 
sure  of  580  psia  (4.00  MPa)  and  a  total  temperature 
of  1350  R  (750  K). 

Experimental  data  includes  surface  pressure  taps 
(24)  and  thermocouples  (32),  mean  profiles  using 
pitot-pressure  and  total  temperature  probes,  and  in¬ 
stability  waves  using  calibrated  constant-current  hot¬ 
wire  anemometry.  Hot-wire  spectra  were  recorded  at 
one-inch  (2.54-cm)  intervals  between  10  and  31  inches 
(25.4  and  78.7  cm),  with  additional  measurements  at 
33,  35,  and  37  inches.  Eour  hot-wires  were  used  for 
these  measurements,  each  calibrated  separately.  There 
is  no  data  for  transition  onset  on  the  blunt  cone,  ac¬ 
cording  to  Stetson. 

AEDC  carried  out  spectral  analyses  of  the  hot-wire 
data  and  printed  out  tables  of  amplitude  vs  frequency 
at  these  streamwise  stations.  These  detailed  exper¬ 
imental  data  are  now  available  on  the  RTO  WG  10 
SG  3  website  (see  above).  Eor  a  detailed  discussion, 
see  Schneider®  and  Kimmel.^® 

Suggested  Approach 

The  following  work  seems  to  be  warranted: 

1.  Use  the  existing  wall-temperature  data  to  perform 
mean-flow  computations  at  isothermal  conditions. 
Since  no  thermocouples  exist  forward  of  7  nose 
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radii,  it  will  be  necessary  to  estimate  the  temper¬ 
ature  distribution  near  the  nose. 

2.  To  obtain  the  experimental  amplification  rates, 
the  experimental  amplitude  data  were  differenti¬ 
ated  by  AEDC.  It  seems  preferable  to  integrate 
the  computational  amplification  rates  to  compare 
amplitude  ratios  vs  frequency  between  compu¬ 
tation  and  experiment.  Tabulated  experimental 
data  is  available  for  this  purpose,  at  the  stream- 
wise  stations  mentioned  above. 

3.  Compare  to  the  extensive  mean-flow  data.  See 
if  there  is  any  consistent  comparison  between  the 
computations  and  experiments. 

4.  Compare  nonlinear  computations  to  the  ampli¬ 
fication  rates  of  the  primary  and  2nd-harmonic 
frequencies  in  the  nonlinear  region. 

Dataset  No.  3:  Sharp  Flared  Cone  at  Mach  6 

Prepared  by 

5.  Schneider 

Introduction 

Three  instability  experiments^^^^^  were  carried  out 
in  the  NASA  Langley  Mach  6  quiet  tunnel  before  it  was 
decommissioned.  Due  to  schedule  limitations,  this  is 
the  only  one  with  calibrated  mean-flow  profiles.  The 
experiments  were  carried  out  on  both  sharp  and  blunt 
cones,  but  extensive  data  is  available  only  for  the  sharp 
cone.  All  measurements  were  made  using  natural  fluc¬ 
tuations,  without  controlled  disturbance  generators. 

The  first  study  that  involved  stability  measurements 
of  a  laminar  boundary  layer  in  a  quiet  hypersonic 
tunnel  was  conducted  by  Lachowicz  and  Chokani  in 
the  early  1990s.  As  the  constant  voltage  anemometer 
(CVA)  was  used  to  operate  the  hot-wire,  and  the  CVA 
was  poorly  understood  the  stability  analysis  was  de¬ 
rived  from  uncalibrated  measurements.  More  recently 
our  knowledge  of  the  CVA  has  significantly  improved; 
fortuitously  also  Lachowicz  and  Chokani  did  also  con¬ 
duct  one  systematic  with  a  calibrated  hot-wire.  In 
the  calibration,  the  hot-wire  was  operated  at  twelve 
overheats  for  seven  mass  fiux/total  temperature  com¬ 
binations.  Then  the  calibrated  hot-wire  was  used  to 
examine  the  disturbance  evolution  on  a  sharp  tipped 
cone  flare  model.  On  the  model  seventeen  streamwise 
stations,  each  with  thirteen  points  across  the  boundary 
layer  were  examined.  At  each  measurement  point  the 
hot-wire  is  operated  at  seven  overheats.  This  paper^* 
describes  and  compares  two  methods  for  obtaining  the 
calibrated  data.  In  addition  the  effect  of  wire  Reynolds 
number,  overheat  ratio,  etc  are  also  discussed.  The 
paper  shall  be  accompanied  by  a  data  report  that  tab¬ 
ulates  the  mean  and  rms  profiles  of  mass  flux  and  total 
temperature.  These  tabulated  data  may  be  useful  for 
purposes  of  code  validation. 


The  experiments^^^^®  focused  on  the  growth  of  first 
and  second-mode  instability  waves  on  an  axisymmet- 
ric  cone  under  cold-flow  conditions,  in  a  quiet  tunnel. 
An  adverse  pressure  gradient  exists  on  the  rear  half  of 
the  model.  Second-mode  amplification  was  observed. 
Gortler  interactions  are  possible  on  the  concave  flare, 
although  no  experimental  evidence  is  available. 

Lachowicz  et  report  agreement  of  2-5%  with 
linear-instabilty  N-factor  computations  for  the  second¬ 
mode  wave  growth  {e.g.,  Fig.  9  in  Lachowicz  et 
However,  the  ratio  of  CVA  fluctuating  voltages  was 
taken  as  the  ratio  of  wave  amplitude.  This  assumption 
was  used  by  Stetson  and  others  for  constant-current 
or  constant-temperature  anemometry,  when  the  mea¬ 
surements  were  carried  out  with  the  hot  wire  at  iden¬ 
tical  mean-flow  conditions. 

If  this  assumption  can  be  verified,  and  any  errors 
caused  by  it  can  be  quantified,  this  dataset  may  be 
the  best  of  the  three.  Although  the  aft  end  of  the 
cone  was  outside  of  the  fully  quiet  region,  this  is  the 
only  data  obtained  under  quiet  conditions.  Static  pres¬ 
sure  data  was  obtained  at  three  azimuthal  positions 
to  check  angle  of  attack,  although  all  boundary-layer 
measurements  were  again  obtained  on  only  one  ray. 
The  mean  flow  profiles  are  in  fairly  good  agreement 
with  existing  computations. 

These  data  were  also  obtained  under  equilibrium- 
wall  temperature  conditions.  That  is,  the  measure¬ 
ments  were  obtained  only  after  the  model-wall  ther¬ 
mocouples  reached  steady-state  temperatures.  This 
occurred  approximately  15-20  minutes  after  initiation 
of  Mach-6  flow,  following  a  subsonic  preheat.  Unfortu¬ 
nately,  as  for  the  Stetson  data,  this  is  not  the  same  as 
adiabatic-wall  conditions,  for  Chokani  states  that  the 
heat-transfer  within  the  model  is  not  negligible. 

Description  of  Experiment 

The  model  was  a  50.8  cm  sharp  cone.  The  first 
25.4  cm  formed  a  straight  cone  with  a  5  deg  half-angle. 
This  was  followed  by  a  gentle  flare  with  a  2.364  m  ra¬ 
dius  of  curvature.  The  profile  shape  is  continuous;  only 
the  second  derivative  is  discontinuous  at  the  match 
point.  The  tip  had  a  nominal  radius  of  2.5  microns. 
The  cone  angle  of  attack  was  nearly  zero;  however, 
most  measurements  were  carried  out  at  a  yaw  angle 
of  approximately  0.1  deg.  and  a  pitch  angle  of  ap¬ 
proximately  0.1  deg  (see  Lachowicz  et  p.  20). 
Appendix  D  of  Lachowicz  et  reports  additional 
measurements  carried  out  after  additional  efforts  were 
made  to  align  the  model  at  zero  angle  of  attack. 

The  measurements  were  carried  out  at  Mach  5.91, 
with  a  total  pressure  of  896  kPa  and  a  total  temper¬ 
ature  of  450K.  The  freestream  Reynolds  number  was 
9.25  X  10®/m  (note  the  typographical  error  in  Lachow¬ 
icz  et  p.  2497). 

Experimental  data  includes  surface  pressure  taps 
(29)  and  thermocouples  (51),  mean  flow  profiles  using 
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hot  wires  and  CVA.  Fluctuating  profiles  were  obtained 
using  same  hot  wires,  but  no  calibrated  fluctuations 
have  so  far  been  reported.  See  Riemann^*  for  the  lat¬ 
est  information. 

Suggested  Approaeh 

It  is  recommended  that  the  reader  contact  Chokani 
or  Schneider  regarding  the  details  of  the  experimental 
conditions  and  the  available  data.  Further  indepen¬ 
dent  computations  are  needed,  to  see  if  the  excellent 
agreement  observed  by  Lachowicz  et  al  can  be  re¬ 
peated  by  others.  Computations  have  been  performed 
by  Balakumar  and  Malik^®  and  Pruett  and  Chang. 
The  computational  approximations  need  to  be  evalu¬ 
ated.  For  example,  computations  should  be  carried 
out  with  an  isothermal  wall-temperature  distribution 
adapted  from  the  experimental  data,  to  determine  the 
effects  of  this  boundary  condition. 

2.  REAL  GAS  FLOWS 

TEAM  LEADER:  G.  CANDLER 

Dataset  No.  1:  Circular  Cylinder 

Prepared  by 

W.  Beck,  G.  Candler 

Introduetion 

Hypersonic  propelled  flight  will  occur  at  Mach  num¬ 
ber  from  about  four  up  to  (perhaps)  fifteen.  Around 
the  body  of  the  vehicle  at  higher  Mach  numbers  shocks 
are  generated,  behind  which  the  air  will  be  dissociated 
to  different  degrees,  depending  not  only  on  the  Mach 
number  but  also  on  the  shock  angle.  Shock-shock 
interactions  {e.g.,  forebody- wing)  can  lead  to  locally 
very  high  heating  rates.  Shock  boundary  layer  inter¬ 
actions,  especially  in  the  regions  of  body  flaps,  lead 
to  influences  on  heating  rates  and  body  flap  efficiency. 
The  physical  state  and  turbulence  level  of  hot  bound¬ 
ary  layers,  especially  on  the  vehicle  forebody  upstream 
of  its  inlet  and  in  its  combustion  chamber,  are  not  fully 
understood.  If,  in  addition  to  the  above,  the  role  of 
high  temperature  effects  (especially  chemical  reaction) 
are  added,  then  the  picture  becomes  immensely  diffi¬ 
cult  to  measure  experimentally  and  to  calculate  using 
CFD.  The  “real”  vehicle  case  (fully  3-D  geometry,  pos¬ 
sibly  transitional  flow  or  turbulence,  full  chemistry, 
complicated  flow  interactions)  is  too  difficult  a  test 
case  for  examining  the  role  of  real  gas  effects. 

Flows  involving  high  temperature  effects,  especially 
kinetics,  are  generally  high  enthalpy,  high  temperature 
flows  which  are  difficult  to  generate  in  a  wind  tunnel, 
which  are  even  harder  to  characterize  using  the  various 
available  classical  and  optical  diagnostic  techniques, 
and  which  are  complicated  and  not  always  amenable 
to  CFD  calculations.  Some  (but  not  all)  of  the  difficul¬ 
ties  encountered  in  carrying  out  aero  (thermo)  dynamic 
testing  and  in  interpreting  the  data  are:  very  short  test 


times,  contaminated  flows  (ablation  material),  nonsta¬ 
tionary  effects  (nozzle  start-up  processes),  driver  gas 
arrival  (in  shock  tunnels),  incomplete  characterization 
of  the  free  stream  (including  knowledge  of  boundary 
layer  transition),  lack  of  accredited  (qualified)  kinetic 
data  for  CFD  calculations,  etc. 

There  are  various  kinetic  schemes  which  have  been 
used  and  further  developed  over  recent  years  to  exam¬ 
ine  air  chemistry  in  high  enthalpy  flows.  The  valida¬ 
tion  of  these  schemes,  with  the  consequent  knowledge 
of  having  the  correct  scheme  at  the  given  conditions 
of  temperature  and  pressure  (be.,  density)  is  consider¬ 
ably  hampered  by  a  lack  of  the  type  of  high  enthalpy 
data  that  is  needed  to  distinguish  between  the  differ¬ 
ent  kinetic  schemes.  For  example,  it  is  easy  to  mea¬ 
sure  stagnation  pressures  on  a  model,  but  this  result 
is  virtually  useless  in  distinguishing  between  kinetic 
schemes.  More  data  are  needed  that  probe  the  flow 
itself,  preferably  at  the  molecular  level. 

The  approach  suggested  here  is  to  simplify  the  prob¬ 
lem  as  much  as  possible  in  order  to  focus  on  the  high 
temperature  effects.  A  circular  cylinder  flow  has  sev¬ 
eral  advantages:  line-of-sight  symmetry  (necessary  for 
optical  techniques),  no  complicated  flow  interactions 
(where  the  physics  are  only  poorly  understood,  e.g., 
unsteadiness,  turbulence,  transition,  etc.),  larger  A/D 
(where  A  is  the  shock  standoff  distance  and  D  is  the 
cylinder  diameter)  than  a  sphere  which  simplifies  the 
use  of  optical  field  measurement  techniques,  and  quite 
straightforward  to  carry  out  CFD. 

Deseription  of  Experiments 

The  experiments  are  in  progress  at  HEG  in  DLR. 
The  test  case  is  a  circular  cylinder  with  diameter  90 
mm  and  length  (including  shaped  holder)  380  mm, 
placed  with  its  axis  transverse  to  the  flow.  The  bow 
shock,  generated  in  front  of  a  circular  cylinder  placed 
in  a  flow  with  high  total  enthalpy,  is  strong  enough 
to  cause  temperatures  in  the  shock  layer  which  are 
sufficient  to  bring  about  dissociation  of  oxygen  and 
molecular  nitrogen.  Depending  on  the  binary  scal¬ 
ing  parameter  pL  (where  p  is  the  density  and  L  is 
the  characteristic  length) ,  the  chemical  relaxation  pro¬ 
cesses  will  either  reach  equilibrium  within  the  shock 
layer  or  will  still  be  reacting  (be.,  remain  in  chemical 
non-equilibrium).  Furthermore,  the  thermal  behav¬ 
ior  of  the  gas  can  also  lead  to  non-equilibrium.  Both 
high  temperature  effects,  thermal  and  chemical  non¬ 
equilibrium,  lead  to  a  modification  of  the  gas  proper¬ 
ties  (temperature,  density,  sound  speed,  etc.)  behind 
the  shock,  which  themselves  lead  to  changes  in  shock 
properties  such  as  shock  standoff  distance  and  shock 
shape. 

Measurements  (see  below)  on  the  cylinder  have  been 
carried  out  in  HEG  at  high  and  intermediate  total  en¬ 
thalpies  (conditions  I  and  HI,  «  22  and  «  13  MJ /kg, 
respectively),  at  low  pressures  (nozzle  reservoir  «  35 
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MPa)  and  with  test  gas  air  and  nitrogen.  Tests  also  at 
high  pressures  (100  MPa)  are  still  under  consideration. 

The  nominal  HEG  freestream  data  are  listed  in  Ta¬ 
ble  4  for  air  at  conditions  I  and  III,  based  on  a  full 
Navier-Stokes  turbulent  non-equilibrium  nozzle  flow 
calculation  (K.  Hannemann)  which  uses  experimen¬ 
tal  input  such  as  shock  speed,  reservoir  pressure  and 
nozzle  wall  heat  transfer  in  order  to  calculate  the  noz¬ 
zle  exit  conditions  (K.  Hannemann,  M.  Schneider)  as 
given  below: 


Table  4  Circular  Cylinder 

Reservoir 

Cond  I 

Cond  III 

Pressure  po  (MPa) 

35 

44 

Temperature  Tg  (K) 

9120 

6980 

Specific  enthalpy  (MJ /kg) 

21.7 

12.5 

Freestream 

Pressure  Poo  (Pa) 

660 

790 

Density  Poo  (kg/m®) 

0.0017 

0.0033 

Temperature  Too  (K) 

1140 

810 

Pitot  pressure  pt^  (kPa) 

55.4 

66.1 

Mach  number  Moo 

8.17 

8.14 

Flow  velocity  Uoo  (m/s) 

5940 

4660 

Freestream  mass  fractions 

[N,] 

0.745 

0.731 

m 

0.047 

0.171 

[NO] 

0.0293 

0.059 

[N] 

4.6  X  10-® 

1.3  X  10-1° 

[0] 

0.178 

0.038 

The  chemical  kinetic  scheme  uses  the  Park  model. 
CFD  results  are  validated  by  comparison  with  mea¬ 
sured  values  of  Po,  Poo,  Poo,Too,  Pt2,  <lt2  and  Uoo- 

The  cylinder  is  equipped  with  pressure  and  temper¬ 
ature  sensors  at  various  positions  (angles)  along  its 
arc  length.  Freestream  static  and  Pitot  (small  sphere) 
pressures  and  stagnation  heat  transfer  (small  sphere) 
are  recorded  each  run  for  calibration,  normalization 
and  statistical  purposes.  Flow  visualization  around 
the  cylinder  is  carried  out  with  the  Schlieren  and  holo¬ 
graphic  interferometric  techniques.  The  use  of  further 
techniques  (spectroscopic  techniques  such  as  PLIF) 
must  await  further  development  and  successful  imple¬ 
mentation. 

Dataset  No.  2:  Shock-Shock  Interaction  on 
Cylindrical  Leading  Edge 

Prepared  by 

M.  Holden 

Introduction 

In  his  early  studies  of  shock-shock  interaction,  Ed- 
ney  speculated  that  regions  of  shock-shock  interaction 
should  be  strongly  influenced  by  real  gas  effects  which 
should  significantly  enhance  heat  transfer  and  pressure 
levels  in  the  interaction  region  by  changing  the  shock 
structure  in  the  shock-shock  interaction.  However, 


preliminary  studies  in  high  enthalpy  flows  at  GAL- 
CIT  and  in  the  HEG  facility  in  Germany  indicated 
that  real  gas  effects  caused  a  decrease  in  both  the  heat 
transfer  and  pressure  amplification  factors.  The  low 
spatial  and  temporal  resolution  of  the  instrumentation 
employed  in  these  studies,  coupled  with  the  data  re¬ 
duction  methods  may  have  had  a  significant  influence 
on  the  accuracy  of  this  data.  Again,  the  question  of  the 
laminar  or  transitional  nature  of  the  phenomena  has 
also  been  raised.  Subsequently,  experimental  studies 
of  shock-shock  interaction  phenomena  were  conducted 
in  the  LENS  facility  employing  models  with  high  spa- 
cial  resolution  and  a  sufficient  frequency  response  to 
accurately  define  the  distribution  of  heat  transfer  and 
pressure  in  the  interaction  region.  Measurements  were 
made  in  both  air  and  nitrogen  flows  and  these  studies 
indicated  that  real  gas  effects  resulted  in  a  signifi¬ 
cant  enhancement  of  the  amplification  factors  in  the 
interaction  regions.  These  studies  were  followed  by 
measurements  on  a  cone-cone  configuration  which  also 
indicated  increased  amplification  factors  as  a  result  of 
real  gas  effects  in  regions  of  shock-shock  interaction. 


Description  of  Experiments 

The  configuration  is  shown  in  Fig.  1.  An  inclined  flat 
plate  generates  a  two  dimensional  oblique  shock  wave 
which  interacts  with  the  bow  shock  generated  by  the 
1.5  in  radius  cylinder.  The  characteristic  of  the  shock- 
shock  interaction  depends  on  a  variety  of  parameters 
including  the  relative  shock  strengths  and  positions  of 
the  oblique  and  blunt  body  shocks  and  these  effects 
were  explored  in  the  experimental  program.  Experi¬ 
ments  were  performed  for  cylinder  radii  of  0.375  and 
1.5  in.  The  experimental  data^^“^®  includes  surface 
temperatures,  heat  transfer  and  pressure  distributions. 
The  freestream  conditions  are  shown  in  Table  5  and 
the  reference  CD-ROM.^® 


Table  5  Shock-Shock  Interaction 


Moo 

Re 

Driven 

Enthalpy 

Run 

(ft-i) 

Gas 

(MJ/kg) 

No. 

8.615 

1.386  X  10® 

Air 

11.3 

50 

9.278 

1.612  X  10® 

Na 

10.5 

48 
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Dataset  No.  3:  Blunt  25°  Cone/60°  Flare 
(Indented  Nose)  Configuration 

Prepared  by 

M.  Holden 

Introduction 

Earlier  experimental  studies  with  indented  nose  tips 
by  Holden  and  more  recently  computations  performed 
by  Candler  have  suggested  that  the  flow  structure  over, 
and  the  aerothermal  loads  generated  on  a  blunted 
cone-cone  configuration,  similar  to  that  of  an  ideal¬ 
ized  indented  nose  tip,  can  be  sensitive  to  real  gas 
effects.  Originally  indented  nose  shapes  were  studied 
because  of  an  interest  in  the  development  of  ablated 
nose  tips,  similar  to  those  generated  with  graphite  nose 
tips  during  ballistic  reentry.  Later  Candler,  in  nu¬ 
merical  studies  designed  to  select  configurations  which 
were  sensitive  to  real  gas  effects,  suggested  a  cone- 
cone  configuration  similar  to  those  studied  earlier  by 
Holden.  To  further  investigate  real  gas  effects  on 
flowfield  structure  and  surface  properties,  Holden  de¬ 
signed  a  blunted  cone/flare  configuration  which  was 
highly  instrumented  with  high-frequency  heat  transfer 
and  pressure  gages.  Measurements  were  attained  to 
provide  distribution  of  the  surface  properties  as  well 
as  information  on  any  high-frequency  flow  oscillations 
that  may  arise.  Test  conditions  weres  selected  to  ob¬ 
tain  air  and  nitrogen  freestream  flows  at  enthalpies  of 
5  and  10  MJ/kg  at  reservoir  pressures  of  270  to  500 
bars.  Computations  employing  the  Navier-Stokes  code 
have  been  made  by  Candler  and  Gnoffo  to  examine 
both  these  flow  configurations.  These  computations 
suggested  that  the  flows  were  unstable,  a  feature  not 
observed  in  the  experiment.  The  numerical  results 
were  also  highly  grid-sensitive  and  caused  significant 
debate  and  controversy  in  the  Navier-Stokes  solver 
community. 

^ -  6.805 - ^ 

1^  4.000^ 


Fig.  2  Blunt  cone  (60°)  (dimensions  in  inches) 


Description  of  Experiments 

The  experiments  were  performed  at  Calspan  -  Uni¬ 
versity  at  Buffalo  Research  Center  (CUBRC)  and  are 
described  in  Holden. The  configuration  is  a  25° 
cone  with  a  60°  flare,  and  is  shown  in  Fig.  2.  The 
experimental  data  includes  surface  temperatures,  heat 
transfer  and  pressure  distributions.  The  freestream 
conditions  are  shown  in  Table  6  and  the  reference  CD- 
ROM.23 

Table  6  Sharp-nose  Cone-Flare 


Afoo 

Re 

Driven 

Enthalpy 

Run 

(ft 

-1 

) 

Gas 

(MJ/kg) 

No. 

8.439 

1.481 

X 

10® 

Air 

11.4 

43 

8.623 

1.379 

X 

10® 

Air 

11.3 

46 

9.279 

1.563 

X 

10® 

N2 

10.5 

45 

10.49 

3.791 

X 

10® 

N2 

5.8 

44 

10.60 

3.774 

X 

10® 

Air 

5.5 

42 

Dataset  No.  4:  Large  Spherically  Blunt-Nose 
Cone- Flare 

Prepared  by 

M.  Holden 

Introduction 

Two  sets  of  studies  were  conducted  with  the  large 
blunt  nose  cone-flare  configuration.  This  model  was 
highly  instrumented  with  heat  transfer  and  pressure 
gages  to  obtain  detailed  surface  distributions  for  a 
range  of  freestream  conditions  at  velocities  between 
10,000  ft/sec  and  15,000  ft/sec  for  both  nitrogen  and 
air  flows.  The  electron  beam  and  the  associated  pneu¬ 
matic  and  electrical  apparatus,  together  with  photo¬ 
multiplier  detection  systems,  were  also  installed  into 
the  model.  A  higher  density  of  surface  instrumenta¬ 
tion  was  employed  in  the  immediate  vicinity  of  the 
cone-flare  junction  to  detect  and  quantify  the  presence 
of  a  laminar  separated  region.  The  flare  region  was 
also  densely  instrumented  to  accurately  quantify  the 
reattachment  compression  region.  Measurements  of 
the  distribution  of  heat  transfer  and  pressure  over  the 
model  were  obtained  for  total  enthalpy  conditions  of  5 
MJ/kg  and  10  MJ/kg,  with  reservoir  pressure  of  500 
bars,  with  both  nitrogen  and  air  used  as  the  test  gas. 
These  conditions  were  selected  to  explore  the  effects  of 
the  Reynolds  number,  which  influences  the  occurrence 
of  transition  and  total  temperature  and  pressure  and 
influences  the  flow  chemistry. 

Description  of  Experiments 

The  configuration  is  a  15°  cone  with  a  30°  flare  and 
an  indented  nose,  which  is  presented  in  Fig.  3.  The 
experimental  data^^’^®  includes  surface  temperatures, 
heat  transfer  and  pressure  distributions.  Flowfield  sur¬ 
veys  were  also  performed  to  examine  real-gas  effects  on 
shock  position  close  to  the  sphere/cone  junction.  The 
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Fig.  3  Large  spherically  blunt  nose  cone-flare 

freestream  conditions  are  shown  in  Table  7  and  the 
reference  CDROM.^® 


Table  7  Blunt-nose  Cone-Flare 


Moo 

Re/ ft 

Driven 

Enthalpy 

Run 

(xlO®) 

Gas 

(MJ/kg) 

No. 

8.508 

1.474 

Air 

11.2 

61 

9.404 

1.519 

N.2 

10.5 

60 

10.55 

3.920 

N2 

5.7 

59 

10.58 

4.083 

Air 

5.4 

63 

3.  LAMINAR  VISCOUS-INVISCID 
INTERACTIONS 
TEAM  LEADER:  M.  HOLDEN 

Dataset  No.  1:  Hollow  Cylinder  Flare 

Prepared  by 
T.  Pot,  B.  Chanetz. 

Introduction 

This  test  case  is  proposed  for  Navier-Stokes  solvers 
validation  in  the  field  of  high  speed  flows  characterized 
by  axisymmetric  strong  shock- wave  boundary-layer  in¬ 
teraction  in  the  fully  laminar  flow  regime.  The  exper¬ 
iment  was  carried  out  in  the  R5Ch  blow-down  Hyper¬ 
sonic  Wind  Tunnel  in  the  Fundamental  and  Experi¬ 
mental  Aerodynamics  Department  of  ONERA  headed 
by  Prof.  J.  Delery,  France.  Further  numerical  stud¬ 
ies  (code-to-code  comparison)  and  comparison  with 
experimental  data  have  been  performed  once  in  the 
frameworks  of  AGARD  Working  Group  18,  then  of 
the  first  Europe-US  data  base  and  the  first  Eastern- 
Western  High  Speed  Flow  Field  data  base.  Details  are 
presented  in  Ghanetz  et 

Description  of  Experiments 
The  model  is  comprised  of  a  hollow  cylinder  with  an 
external  diameter  of  65  mm  and  an  interned  diameter 
of  45  mm  as  shown  in  Fig.  4.  The  leading  edge  is  sharp 
with  an  angle  of  15°.  A  30°  flare  is  located  101.7  mm 


downstream  of  the  cylinder  leading  edge.  The  flare 
is  followed  by  a  cylindrical  part  having  a  diameter  of 
115  mm  and  a  length  of  25  mm.  The  advantage  of 
this  configuration  is  twofold:  i)  to  put  baseflow  far 
enough  downstream  of  the  interaction  region,  and  ii) 
to  eliminate  three  dimensional  effects.  The  origin  of 
the  coordinate  system  is  located  at  the  leading  edge  of 
the  hollow  cylinder.  The  mesh  coordinates  are  nondi- 
mensionalized  by  the  reference  length  L  =  0.1017  m. 
The  test  conditions  (Table  8)  are  perfect  gas  (7  =  1.4), 
axisymmetric,  fully  laminar  flow,  laminar  viscosity  by 
Sutherland’s  Law.  Experimental  data  includes  surface 
pressure  coefficient  distribution,  surface  heat  transfer 
(Stanton  number)  distribution,  surface  oil-flow  visu¬ 
alization  (skin  friction  lines) ,  fiowfield  visualization  by 
the  Electron  Beam  Fluorescence  (EBF)  technique,  and 
density  profiles. 


Fig.  4  Hollow  cylinder  flare 

Table  8  Hollow  Cylinder  Flare 

Quantity  Value 

Moo  9.91 

Re/m  1.86  x  10® 

Too  51  K 

Poo  6.3  Pa 

^wall  293  K 

L  0.1017  m 


Dataset  No.  2:  Hollow  Cylinder  Flare 

Prepared  by 

M.  S.  Holden 

Introduction 

This  study  was  conducted  to  obtain  detailed  exper¬ 
imental  measurements  of  the  heat  transfer  and  pres¬ 
sure  distribution  over  two  cylinder-flare  models,  one 
of  which  replicates  a  configuration  employed  earlier  in 
French  code  validation  experiments  (see  above)  and 
a  second  which  employs  a  longer  flare  to  obtain  a 
better  defined  downstream  boundary  condition.  The 
boundary  layer  was  laminar  over  the  entire  cylinder- 
flare  surface  in  all  cases,  thereby  enabling  a  detailed 
comparison  between  DSMG,  Navier  Stokes  and  hybrid 
codes.  A  detailed  description  of  the  experiments  is 
presented  in  Holden  and  Wadhams.^® 
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Fig.  5  Hollow  cylinder  with  short  flare  (dimen¬ 
sions  in  inches) 


Fig.  6  Hollow  cylinder  with  long  flare  (dimensions 
in  inches) 

Description  of  Experiments 

The  two  model  configurations  are  shown  in  Figs.  5 
and  6.  The  first  configuration  employs  a  shorter  fiare 
configuration  which  is  identical  in  geometry  to  the 
configuration  used  in  earlier  French  code  validation 
experiments.  Because  we  had  believed  that  the  sep¬ 
arated  region  on  this  configuration  is  infiuenced  by 
the  expansion  at  the  downstream  end  of  the  conical 
fiare,  we  had  requested  P.  Gnoffo  (NASA  Langley) 
perform  calculations  to  determine  the  length  of  fiare 
required  to  ensure  that  the  separated  region  is  not  ef¬ 
fected  by  the  shoulder  and  region  of  constant  pressure 
generated  upstream  of  the  end  of  the  fiare.  Based  on 
results  of  Gnoffo’s  calculations,  we  selected  the  sec¬ 
ond  of  the  two  configurations  (Fig.  6),  which  has  a 
fiare  of  sufficient  length  so  that  the  separated  region 
is  uninfluenced  by  the  downstream  end  of  the  conical 
fiare.  These  two  model  configurations  were  employed 
in  studies  at  Calspan  at  the  test  conditions  indicated  in 
Table  9.  The  test  conditions  were  selected  on  the  basis 
of  many  fiowfield  calibration  studies  to  ensure  stream- 
wise  and  lateral  fiowfield  uniformity.  The  experimental 
measurements  include  detailed  high  frequency  surface 
temperatures,  heat  transfer  and  pressure  distributions 
and  Schlieren  photographs  of  the  flow. 

The  fiowfield  physics  for  this  study  consist  of  perfect 
gas,  axisymmetric  fiowfield,  fully  laminar  shock  wave¬ 
boundary  layer  interactions,  and  flow  separations.  The 
test  gas  is  nitrogen. 


Table  9 

Hollow  Cylinder 

Flare 

Case 

Moo 

Re/ ft 
(xlO^) 

Short 

1 

9.30 

10.8 

2 

9.44 

8.33 

3 

11.3 

6.28 

4 

11.4 

10.7 

5 

11.5 

7.92 

Extended 

1 

9.33 

10.9 

2 

9.38 

6.51 

3 

9.45 

7.89 

4 

9.50 

5.60 

5 

11.1 

4.75 

6 

11.4 

7.19 

7 

11.4 

7.89 

8 

11.4 

10.7 

Dataset  No.  3:  Sharp-Nose  and  Blunt-Nose 
Cone-Flare  (Indented  Nose)  Configuration 

Prepared  by 

M.  Holden,  G.  Candler,  J.  Harvey 
Introduction 

The  flows  around  hypersonic  vehicles  will  have  re¬ 
gions  of  shock/shock  and  shock/boundary  layer  inter¬ 
action  which  can  induce  embedded  regions  of  sepa¬ 
rated  flows  over,  for  example,  control  surfaces,  and 
in  and  ahead  of  airbreathing  ramjet  or  scramjet  en¬ 
gines.  Oblique  shock  waves  generated  inside  the  engine 
can  also  cause  extensive  regions  of  three  dimensional 
boundary  layer  separation. 

Recent  work  has  shown  that  the  flow  generated  by 
a  cone-flare  geometry  is  very  challenging  to  compute. 
This  flow,  with  an  attached  leading  edge  shock  and  a 
detached  bow  shock  from  the  second  cone,  produces  a 
complicated  interaction  between  the  shock  waves  and 
the  separation  zone.  The  size  of  the  separated  region 
is  very  sensitive  to  the  shock  angles  and  the  strength  of 
the  shock  interaction.  Because  the  size  of  separated  re¬ 
gion  can  be  measured  with  surface  measurements  and 
with  flow  visualization  methods,  this  fiowfield  is  an 
excellent  test  case  for  code  validation.  Also,  because 
the  flow  is  axisymmetric,  the  flowfields  require  only 
two-dimensional  simulations. 

The  cone-flare  configuration  is  a  challenge  for  CFD 
because  it  is  critical  to  resolve  the  initial  boundary 
layer  growth  near  the  model  leading  edge.  Addition¬ 
ally,  the  shock  interaction  must  be  captured  accurately 
in  order  to  correctly  predict  the  interaction  between 
the  shock  waves  and  the  separation  zone.  However, 
Wright  et  have  shown  that  sufficiently  fine  grids 
and  careful  grid  generation  CFD  methods  can  repro¬ 
duce  this  type  of  fiowfield.  The  more  challenging 
higher  cone  angle  cases  given  in  Wright  et  al  have  a 
transitional  shear  layer,  and  as  a  result  are  not  suit¬ 
able  for  benchmark  CFD  validation  studies.  Therefore 
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a  new  set  of  experiments  with  purely  laminar  flow 
has  been  designed  and  performed  at  Calspan  by  M. 
Holden.  The  models  are  heavily  instrumented  to  accu¬ 
rately  locate  the  separation  and  reattachment  points. 
The  experiments  were  selected  on  the  basis  of  ear¬ 
lier  studies  in  separated  flows  to  ensure  that  the  flows 
over  the  cone-flare  configurations  will  remain  laminar 
throughout  the  interaction  region. 

h -  7.625 - H 


H —  4.000-H 


Fig.  7  Sharp  cone  (55°)  (dimensions  in  inches) 


7.625 - * 

4.000-^ 


Fig.  8  Sharp  cone  (60°)  (dimensions  in  inches) 

Description  of  Experiments 

Several  of  the  model  configurations  are  shown  in 
Figs.  7,  8,  9  and  10.  A  full  description  of  the  models 
is  provided  in  Holden  and  Wadhams.^®  The  blunted 
25°  cone  with  a  60°  flare  was  employed  in  earlier  ex¬ 
perimental  studies  with  the  objective  of  investigating 
the  effects  of  real  gas  on  flowfield  geometry.  These 
experimental  studies  indicated  that  the  flowfield  was 
steady,  and  that  there  were  significant  differences  be- 


6.805 - * 

4.000^ 


Fig.  9  Blunt  cone  (55°)  (dimensions  in  inches) 


6.805 - - 

1^  4.000^ 


Fig.  10  Blunt  cone  (60°)  (dimensions  in  inches) 

tween  the  heat  transfer  and  pressure  distributions  for 
nitrogen  and  air  flows.  However,  computations  by  G. 
Candler  and  P.  Gnoffo  indicated  that  the  flows  over 
this  configuration  are  unsteady,  and  grid  refinement 
studies  indicated  the  solutions  were  grid  dependent. 
We  repeated  the  experiment  at  much  lower  Reynolds 
numbers  (see  Table  10)  to  ensure  that  the  flow  remains 
completely  laminar  throughout  the  interaction  region 
over  this  configuration. 

We  investigated  the  effect  of  nose  bluntness  by  em¬ 
ploying  a  sharp  25°  nosetip.  In  addition,  we  performed 
studies  with  the  blunt  and  sharp  nosetips  with  a  55° 
conical  flare  in  place  of  the  original  60°  flare.  The 
55°  flare  configuration  was  selected  in  conjunction  with 
calculations  done  by  P.  Gnoffo,  which  indicated  that 
the  basic  flow  over  the  configuration  was  stable,  al¬ 
though  there  appeared  to  be  some  problems  in  regards 
to  instability  in  the  recirculation  region. 

A  bow  shock  is  generated  over  the  sharp  or  spher- 
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ically  blunted  nosetip  which  interacts  with  the  de¬ 
tached  bow  shock  generated  over  the  flare.  A  region 
of  shock/shock  interaction  is  generated  at  the  inter¬ 
section  of  the  two  shocks  and  the  resultant  shock/- 
boundary  interaction  induces  an  embedded  separated 
region  at  the  cone/flare  junction.  The  coupling  be¬ 
tween  the  shock/shock  interaction  and  the  separated 
zone  makes  this  flow  very  sensitive  to  grid  refinement 
in  numerical  simulations.  The  computed  flowfleld  is 
also  very  sensitive  to  the  grid  spacing  at  the  nosetip 
and  in  the  reattachment  compression  region  at  the 
base  of  the  shock/shock  interaction. 

Table  10  Sharp-  and  Blunt-Nose  Cone-Flare 


Cone-Flare 

Nose 

Moo 

Re/ ft 

Radius  (in) 

(xlO^) 

25°/55° 

sharp 

9.42 

6.18 

9.48 

7.79 

9.59 

4.25 

11.3 

4.76 

0.250 

9.42 

6.27 

9.43 

7.26 

9.56 

4.34 

11.3 

4.33 

0.288 

9.42 

6.41 

11.5 

4.18 

0.600 

9.47 

3.80 

25°/60° 

sharp 

11.4 

4.70 

0.250 

11.3 

4.44 

The  experimental  data  include  detailed  heat  transfer 
and  pressure  distributions.  The  55°  model  is  instru¬ 
mented  with  20  pressure  gauges  and  34  temperature 
gauges  while  the  60°  model  is  instrumented  with  18 
pressure  gauges  and  33  temperature  gauges  in  the  sep¬ 
aration  and  reattachment  locations.  Measurements 
were  performed  with  both  the  blunt  and  sharp  double 
cones  for  flare  angles  of  55°  and  60°.  The  freestream 
conditions  are  shown  in  Table  10.  The  test  gas  is  ni¬ 
trogen. 

4.  SHOCK-SHOCK  INTERACTIONS 
TEAM  LEADERS:  J.  SCHMISSEUR 
S.  WALKER 


and  Experimental  Aerodynamics  Department  of  ON- 
ERA  headed  by  Prof.  J.  Delery,  France.  Details  are 
presented  in  Pot  et  Moss  et  Delery  and 
Chanetz®®  and  Reijasse  et 

Description  of  Experiments 

The  arrangement  of  the  cylinder  and  the  wedge 
models  in  the  test  section  is  shown  in  Fig.  11.  The 
relative  locations  of  the  cylinder  and  the  wedge  are 
selected  in  order  to  obtain  a  Type  IV  shock-shock  in¬ 
teraction.  The  wedge  angle  is  20°  and  the  wedge  length 
is  equal  to  Lwedge  =  0.1  m.  Its  span  is  0.1  m.  The 
axis  of  the  cylinder  is  perpendicular  to  the  free  stream 
direction.  The  distance  between  the  leading  edge  of 
the  shock  generator  and  the  cylinder  center  is  0.11  m 
along  the  x-axis  (longitudinal  parallel  to  the  upstream 
flow  direction)  and  0.053  m  along  y-ajds  (perpendicu¬ 
lar  to  the  upstream  flow  direction).  The  origin  of  the 
coordinate  system  is  located  at  the  nose  of  the  cylin¬ 
der.  The  mesh  coordinates  are  nondimensionalized  by 
the  cylinder  radius  =  0.008  m. 


Table  11  Shock-Shock  Interaction 

Quantity  Value 

Moo  9.91 

Re/m  1.86  x  10® 

Too  51  K 

Poo  6.3  Pa 

^wall  293  K 

L  0.1017  m 


Dataset  No.  1:  Shock-Shock  Interaction 

Prepared  by 
T.  Pot,  B.  Chanetz 

Introduction 

This  problem  has  been  proposed  to  evaluate  CFD 
codes  for  the  analysis  of  complex  shock-shock  inter¬ 
actions.  The  interaction  between  an  impinging  shock 
wave  generated  by  a  wedge  and  the  bow  shock  around 
a  circular  cylinder  (Type  IV  shock-shock  interaction  in 
the  classification  of  Edney)  is  considered.  The  experi¬ 
ments  were  conducted  in  the  R5Ch  in  the  Fundamental 


The  flowfleld  is  described  as  a  perfect  gas  (j  =  1.4), 
two-dimensional,  fully  laminar  flow,  and  laminar  vis¬ 
cosity  by  Sutherland’s  Law. 

The  cylinder  is  equipped  with  thermocouples  and 
pressure  taps.  The  experimental  data  includes  wall 
pressure,  wall  heat  flux,  flowfleld  visualization  by  the 
Electron  Beam  Fluorescence  technique,  and  density 
and  temperature  measurement  in  the  flowfleld  inter¬ 
action  (more  than  100  data  points). 

The  selected  flow  conditions  are  listed  in  Table  11. 
The  effective  stagnation  conditions  (pt^  =  2.5  x  10® 
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Pa  and  Tt^  =  1050  K)  yield  the  freestream  properties 
shown  in  the  Table. 

Dataset  No.  2:  Shock-Shock  Interactions  For 
Ideal  Gas  Flows 

Prepared  by 
M.  S.  Holden 

Intro  duction 

An  extensive  series  of  studies^®’ was  con¬ 
ducted  over  a  range  of  Mach  numbers  from  10  to  16 
to  define  the  aerothermal  loads  generated  in  regions  of 
shock-shock  interaction  from  the  rarefied  flow  to  the 
fully  continuum  turbulent  flow  regimes.  Detailed  heat 
transfer  and  pressure  measurements  were  made  in  the 
48-inch,  96-inch  and  LENS  shock  tunnels.  The  results 
of  these  studies  were  analyzed  to  provide  guidance  to 
predict  the  heating  enhancement  factors  in  laminar 
transitional  and  turbulent  flow  regimes.  The  test  con¬ 
ditions  presented  here  (Run  Nos.  36,  38,  39,  and  40) 
is  a  subset  of  the  experiments  for  fully  laminar  flows; 
for  all  but  one  test  condition,  where  the  flows  in  the 
interaction  region  are  “fully”  turbulent. 


Description  of  Experiments 

The  model  configuration  is  shown  in  Fig.  12.  The 
regions  of  shock-shock  interaction  studied  were  gen¬ 
erated  over  a  series  of  cylindrical  leading  edge  con¬ 
figurations  with  nose  radii  of  0.138,  0.375,  and  1.5 
inches.  Each  of  these  leading  edges  was  densely  in¬ 
strumented  with  heat  transfer  instrumentation  placed 
to  have  a  circumferential  resolution  within  one  de¬ 
gree.  The  thin- film  instrumentation  was  deposited  on 
a  low-conductivity  surface  to  minimize  measurement 
errors  associated  with  lateral  conduction  in  the  large 
heat  transfer  gradients  generated  in  the  region  of  peak 
heating.  The  high-frequency  response  of  the  thin-fllm 
instrumentation  was  also  a  key  factor  in  accurately 
determining  the  heating  distribution  for  shock-shock 
interactions,  which  exhibited  intrinsic  flow  unsteadi¬ 
ness. 

The  flowfleld  physics  for  this  study  consist  of  perfect 
gas,  planar  flowfleld,  shock-shock  interactions,  fully 
laminar  and  highly  turbulent  flows. 

The  experimental  data  includes  surface  tempera¬ 
ture,  heat  transfer  and  pressure  distributions  and 
Schlieren  photographs.  A  sample  of  freestream  test 


conditions  is  shown  in  Table  12.  Additional  informa¬ 
tion  is  available  in  Holden.^® 


Table  12  Shock-Shock  Interaction 


Runs 

36 

38 

39 

40 

Mi 

5.675 

5.534 

5.490 

5.517 

po  X  10“^  (psia) 

4.568 

4.659 

4.097 

4.201 

Ho  X  10“’’  (ftVs^) 

4.686 

4.451 

4.401 

4.439 

To  X  10“®  (°R) 

5.987 

5.784 

5.731 

5.767 

M 

13.99 

14.15 

14.17 

14.15 

U  X  10“®  (ft /sec) 

9.561 

9.319 

9.268 

9.308 

T  (°R) 

193.5 

179.9 

177.2 

179.4 

Poo  X  10^  (psia) 

5.324 

5.312 

4.657 

4.788 

poo  X  lO’^  (slug/ft®) 

2.298 

2.467 

2.195 

2.230 

Poo  X  10^  (slug/ft-sec) 

1.560 

1.449 

1.427 

1.444 

Re  X  10“^  (ff^) 

1.409 

1.587 

1.426 

1.437 

Po  X  10^  (pitot)  (psia) 

1.339 

1.366 

1.202 

1.232 

Dataset  No.  3:  Steady  Shock  Wave  Reflection 
Transition  and  Hysteresis 

Prepared  by 

M.  Ivanov,  A.  Kudryavtsev,  D.  Khotyanovsky 
Introduction 

The  intersection  of  two  symmetric  oblique  planar 
shock  waves  results  in  either  a  Regular  Reflection  (RR) 
or  Mach  Reflection  (MR)  depending  on  the  upstream 
Mach  number  and  shock  angle  as  illustrated  in  Fig.  13. 
At  sufficiently  high  Mach  numbers,  three  regimes  ex¬ 
ist  delineated  by  cejv  (the  von  Neumann  criterion)  and 
an  (the  detachment  criterion)  as  shown  in  Fig.  14. 
For  incident  shock  angle  a  <  Un,  only  RR  is  pos¬ 
sible,  and  for  a  >  an  only  MR  is  possible.  In  the 
region  un  <  Oi  <  an  (the  dual  solution  domain)  both 
types  of  reflection  are  theoretically  possible.  It  was  hy¬ 
pothesized  by  Hornung  that  a  hysteresis  phenomenon 
should  be  observed  when  changing  the  incident  shock 
wave  angle  continuously,  i.e.,  the  transition  from  RR 
to  MR  should  occur  at  a  ^  an  and  reverse  transi¬ 
tion  from  MR  to  RR  at  a  «  ajy.  The  hysteresis  was 
first  obtained  numerically  by  Ivanov  et  al.^^  Recently, 
experiments  by  Ivanov  et  al^^  have  confirmed  the  hy¬ 
pothesis  of  Hornung. 


Regular  Reflection  Mach  Reflection 

Fig.  13  Regular  and  Mach  reflections 

The  test  case  proposed  is  to  perform  unsteady  2D 
Euler  simulations  of  shock  wave  reflection  in  the  dual 
solution  domain,  to  determine  the  angles  of  the  tran- 
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sition  from  RR  to  MR,  and  back,  and  to  obtain  nu¬ 
merically  the  hysteresis  in  simulations  with  increasing 
and  decreasing  the  incident  shock  wave  angle.  At  first 
stage,  only  2D  computations  are  proposed. 

This  problem  is  very  attractive  in  terms  of  verifica¬ 
tion  of  the  available  numerical  codes  for  solving  Euler 
equations.  The  numerical  simulation  of  the  transition 
from  RR  to  MR  is  a  challenging  task  for  various  CFD 
methods  because  of  the  following  features  present  in 
the  flow:  strong  shock  waves,  low  density  regions  be¬ 
hind  the  wedge,  possible  unsteady  phenomena  such  as 
the  instability  of  a  slip  surface  at  MR;  importance  of 
accurate  prediction  of  the  angles  of  the  transition  be¬ 
tween  RR  and  MR;  existence  of  two  stable  states  (RR 
and  MR),  and  hysteresis. 

Description  of  Experiments 

Two  symmetrically-spaced  15°  wedges  are  installed 
in  a  supersonic  flow  at  the  angle  of  attack  9.  Incident 
shock  waves  generated  by  either  wedge  interact  near 
the  plane  of  symmetry.  This  symmetrical  arrangement 
allows  to  eliminate  viscous  effects  inevitable  at  shock 
wave  reflection  from  the  wall.  Figs.  15  and  16  give 
the  schematic  of  the  flow  for  RR  and  MR:  IS  (incident 
shock),  RS  (reflected  shock),  EF  (expansion  fan),  MS 
(Mach  stem),  SS  (slip  surface),  P  (reflection  point), 
and  T  (triple  point).  The  flowfield  is  a  perfect  gas 
(7  =  1.4),  inviscid  two-dimensional  flow,  with  shock 
wave  reflection  transition. 


All  the  lengths  are  non-dimensionalized  by  the 
wedge  length  w.  The  distance  between  trailing  edge 
of  the  wedge  and  the  plane  of  symmetry  is  chosen  to 


be  g/w  =  0.42  to  provide  clear  conditions  when  the 
reflected  shock  at  MR  does  not  impinge  on  the  wedge 
surface  (which  results  in  a  choked  flow)  even  at  large 
incident  shock  angles,  and,  on  the  other  hand,  the  ex- 
pemsion  fan  does  not  interact  with  the  incident  shock. 
The  wedge  is  rotated  around  its  trailing  edge  {i.e.,  g/w 
is  kept  constant),  thus  changing  its  angle  of  attack  6, 
and  correspondingly  the  incident  shock  angle  a.  The 
angles  of  the  transition  uur^mr,  ocmr^rr  between 
RR  and  MR  may  be  compared  with  theoretical  crite¬ 
ria,  and  the  results  computed  by  various  codes.  The 
Mach  stem  heights  at  MR  at  various  a  may  also  be 
compcired  with  experiment. 

Experiments®®  were  conducted  at  ITAM  T313  facil¬ 
ity  at  Moo  =  4  and  5.  At  =  5,  the  corresponding 
theoretical  criteria  are  ajv  =  30.79°,  an  =  39.33°. 
The  experimental  data  include  the  Mach  stem  heights 
measured  at  different  a,  Schlieren  visualizations,  and 
laser  sheet  patterns.  The  relative  spanwise  extent  of 
the  wedge  was  varied  from  b/w  =  0.66  to  b/w  =  3.75 
in  these  experiments.  The  experiments  were  made  at 
the  Reynolds  number  based  on  the  free-stream  param¬ 
eters  and  the  length  of  the  wedge  Re,„  =  2  x  10®,  i.e. 
the  wedge  boundary  layer  effects  can  be  neglected  and 
Euler  simulations  are  justified  in  this  case. 

Previous  simulations  by  Ivanov  et  were  made 

in  the  following  manner.  For  the  incident  shock  wave 
angle  a  nem  emd  slightly  below  the  von  Neumann  con¬ 
dition  Qiv,  a  RR  was  obtained.  Using  this  RR  field 
as  the  initial  data,  the  angle  of  the  wedge  9  was  in¬ 
creased  until  the  transition  from  RR  to  MR  occurred 
at  the  incident  shock  angle  auR^MR  (which  is  near 
and  slightly  above  ou).  After  that,  with  the  ob¬ 
tained  MR  configuration  as  the  initial  data,  the  angle 
of  the  wedge  was  decreased  down  to  ajv,  near  which 
at  a  =  auR^RR  the  transition  to  RR  should  oc¬ 
cur.  The  computations  thus  covered  the  whole  dual 
solution  domain  and  the  hysteresis  loop  was  obtained 
during  decreasing  and  increasing  the  incident  shock 
angles.  The  transition  angles  in  the  computations  are 
very  grid-sensitive,  and  are  also  dependent  on  the  nu¬ 
merical  viscosity.  Mach  stem  heights  correspond  well 
with  the  experimental  ones,  except  at  large  a  where  3D 
effects  become  more  pronounced  in  the  experiments. 
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Dataset  No.  4:  3D  Steady  Shock  Wave  Reflection 

Prepared  by 

M.  Ivanov,  A.  Kudryavtsev,  D.  Khotyanovsky 
Introduction 

The  investigation  of  the  nature  of  the  transition  be¬ 
tween  RR  and  MR  of  steady  shock  waves  between  two 
wedges  is  incomplete  without  a  consideration  of  three- 
dimensional  effects  on  shock  wave  configuration.  A 
finite  spanwise  extension  of  a  test  model  may  have  an 
important  influence  on  the  experimental  results.  Inci¬ 
dent  shock  waves  become  non-planar  because  of  their 
interaction  with  an  expansion  flow  generated  by  lat¬ 
eral  edges  of  the  wedges.  It  is  clear  that,  depending 
on  the  geometrical  sizes,  the  shock  wave  interaction 
may  be  purely  3D  or  contain  a  portion  where  the  inci¬ 
dent  shock  waves  remain  planar.  In  the  case  of  MR,  a 
subsonic  region  behind  the  Mach  stem  exists,  and  as  a 
result,  the  central  plane  flow  should  be  influenced  by 
the  3D  effects.  Strictly  speaking,  the  existing  theoret¬ 
ical  criteria  of  transition  between  RR  and  MR  cannot 
be  applicable  to  real  3D  flows. 

The  conventional  Schlieren  experimental  technique 
gives  flow  images  integrated  along  the  spanwise  direc¬ 
tion  and,  consequently,  does  not  allow  us  to  investigate 
the  3D  structure  of  shock  wave  configurations.  The 
laser  sheet  method  has  been  utilized  by  Ivanov  et 
to  obtain  the  series  of  spanwise  slices  and  reconstruct 
the  whole  3D  flow  patterns. 

Numerical  modelling  can  provide  much  more  de¬ 
tailed  information  concerning  the  3D  structure  of 
shock  wave  reflection.  The  test  case  proposed  is  to 
perform  the  3D  Euler  simulations  of  the  RR  and  MR 
of  steady  shock  waves  generated  by  two  symmetrical 
wedges  of  finite  spanwise  extension.  The  comparison  of 
numerical  and  experimental  data  gives  an  opportunity 
to  validate  numerical  codes  for  simulation  of  super¬ 
sonic  flows  with  complicated  3D  shock  wave  interac¬ 
tions  and  investigate  this  challenging  problem  having 
great  theoretical  and  practical  importance. 

Description  of  Experiments 

Two  symmetrically  spaced  wedges,  which  have  the 
angle  15°,  are  installed  in  a  uniform  supersonic  stream 
at  the  angle  of  attack  9  and  generate  two  interact¬ 
ing  oblique  shock  waves.  This  symmetrical  arrange¬ 
ment  allows  us  to  eliminate  viscous  effects  inevitable 
at  shock  wave  reflection  from  the  wall.  The  flowfleld  is 
a  perfect  gas  (7  =  1.4),  inviscid  3D  flow,  steady  shock 
wave  reflection. 

Fig.  17  gives  the  schematic  of  the  test  model  and 
the  sketch  of  shock  wave  configuration  in  the  vertical 
plane  of  symmetry.  Here  IS  is  the  incident  shock,  RS 
is  the  reflected  shock,  MS  is  the  Mach  stem,  SS  is 
the  slip  surface  emanating  from  the  triple  point  T, 
h  and  g  denote  distances  from  the  horizontal  plane 
of  symmetry  to  the  leading  and  trailing  edges  of  the 


wedge,  respectively,  w  is  the  wedge  length,  b  is  its  span, 
and  s  is  the  half-height  of  the  Mach  stem. 


Fig.  17  3-D  shock  wave  reflection 

Experiments^^’ were  conducted  at  ITAM  T313 
blow-down  wind  tunnel  at  a  freestream  Mach  number 
Moo  =  4-  The  experimental  data  include  the  Mach 
stem  heights  measured  at  different  a,  Schlieren  visu¬ 
alizations,  and  laser  sheet  patterns  taken  in  different 
spanwise  cross-sections. 

The  relative  span  of  the  wedge  in  these  experiments 
was  varied  from  b/w  =  0.66  when  the  incident  shocks 
in  the  vertical  plane  of  symmetry  are  affected  by  the 
lateral  expansions  to  b/w  =  3.75  when  the  shocks  in 
the  central  part  of  the  flow  were  planar.  The  exper¬ 
iments  were  made  at  the  Reynolds  number  based  on 
the  free-stream  parameters  and  the  length  of  the  wedge 
Re^  ~  2  X  10®,  i.e.,  the  wedge  boundary  layer  effects 
can  be  neglected  and  Euler  simulations  are  justified  in 
this  case. 

Previous  Simulations 

Previous  3D  simulations  by  Ivanov  et  al'^^  and 
Kudryavtsev  et  were  performed  both  with  Euler 
code  and  DSMC  method  for  different  flow  parame¬ 
ters  and  geometries.  The  computations  revealed  the 
hysteresis  phenomenon  at  the  transition  from  RR  to 
MR  and  back,  similar  to  that  observed  earlier  in  2D 
case.  The  Mach  stem  heights  in  experiments  and 
computations  are  in  excellent  agreement,  and  strong 
dependence  of  the  Mach  stem  height  on  the  wedge 
span  is  confirmed.  A  close  resemblance  of  numerical 
and  experimental  flow  patterns  was  shown  with  the 
help  of  numerical  Schlieren  pictures  integrated  along 
the  spanwise  coordinate  as  well  as  by  comparison  of 
the  numerical  and  experimental  (laser-sheet)  slices  at 
different  spanwise  locations.  New  and  surprising  de¬ 
tails  of  the  3D  shock  wave  configurations  were  revealed 
at  first  numerically  and  then  confirmed  experimen¬ 
tally.  These  details  include  the  peripheral  MR  when 
the  overall  configuration  is  regular,  the  non-monotonic 
variation  of  the  Mach  stem  height  in  the  spanwise  di¬ 
rection  at  the  MR,  and  the  shock  wave  configuration 
with  intermittent  type  of  reflection:  MR/RR/MR. 


5.  SHOCK  WAVE  TURBULENT 

BOUNDARY  LAYER  INTERACTIONS 

TEAM  LEADER:  D.  KNIGHT 

Prepared  by 

D.  Knight,  A.  Panaras,  A.  Zheltovodov 
Introduction 

The  Conventional  Reynolds-averaged  Navier-Stokes 
(CRANS)  methods"  have  typically  been  unsuccess¬ 
ful  in  predicting  several  of  the  important  aerody¬ 
namic  and  aerothermodynamic  properties  of  shock 
wave  turbulent  boundary  layer  interactions  (also  de¬ 
noted  as  “interactions”)  with  strong  flow  separation. 
CRANS  methods  do  not  accurately  predict  the  ob¬ 
served  mean  surface  pressure,  skin  friction  and  heat 
transfer  in  nominally  two  dimensional  separated  in¬ 
teractions.^®^*  CRANS  methods  do  provide  accurate 
predictions  of  surface  pressure  for  many  moderately 
separated  three  dimensional  interactions,  but  do  not 
provide  an  accurate  prediction  of  skin  friction  and  heat 
transfer  for  separated  interactions. Additionally, 
CRANS  does  not  provide  rms  fluctuating  surface  pres¬ 
sure. 

The  inability  to  accurately  predict  the  desired  aero¬ 
dynamic  and  aerothermodynamic  properties  of  shock 
wave-turbulent  boundary  layer  interactions  signifi¬ 
cantly  degrades  the  capability  for  design  of  effective 
hypersonic  vehicles.  Research  in  the  development  of 
new  methodologies  for  accurate  prediction  of  these 
flows  is  critically  important  to  the  development  of  the 
next  generation  of  hypersonic  vehicle  designs. 

Three  approaches  have  recently  been  identified  as 
promising  candidates  for  improving  prediction  of  shock 
wave  turbulent  boundary  layer  interactions  -  specif¬ 
ically,  for  achieving  improved  prediction  of  aerody¬ 
namic  and/or  aerothermodynamic  properties  com¬ 
pared  to  CRANS  methods.  The  first  approach  is 
Direct  Numerical  Simulation  (DNS)  wherein  the  en¬ 
tire  spectrum  of  turbulent  fluctuations  is  resolved  in 
the  simulation.  The  second  approach  is  Large  Eddy 
Simulation  (LES)  wherein  the  large  energy-containing 
eddies  are  explicity  computed  and  the  small  energy- 
dissipating  scales  modelled.  Promising  results  have  re¬ 
cently  been  obtained  by  Hunt  and  Nixon®°  and  Urbin 
et  for  nominally  two-dimensional  shock  wave  tur¬ 
bulent  boundary  layer  interactions.  The  third  ap¬ 
proach  is  Engineered  Reynolds-averaged  Navier-Stokes 
(ERANS)  wherein  a  CRANS  model  is  modified  to  im¬ 
prove  accuracy  for  a  specific  class  of  flows.  ERANS 
models  may  include  modifications  based  on  knowledge 
of  the  flowfleld  of  shock  wave  turbulent  boundary  layer 
interactions  {e.g.,  Coakley  and  Huang, Panaras®*) 
or  ad  hoc  modifications  for  a  given  configuration  {e.g., 
Bedarev  et  al^"^). 

“The  term  “conventional”  implies  that  there  is  no  problem- 
specific  modification  to  the  general  formulation  of  the  particular 
RANS  model. 


The  objective  of  this  topical  group  is  to  assess  the  ca¬ 
pability  of  DNS,  LES  and  ERANS  to  predict  the  aero¬ 
dynamic  and  aerothermodynamic  properties  of  con¬ 
figurations  with  shock  wave  turbulent  boundary  layer 
interactions  including  strong  separation. 

Dataset  No.  1:  2-D  Supersonic  Compression 
Corner 

Prepared  by 
A.  Zheltovodov 

Description  of  Experiments 

A  two  dimensional  supersonic  equilibrium  turbulent 
boundary  layer  flows  past  a  compression  corner  of  an¬ 
gle  a  as  indicated  in  Eig.  18.  A  shock  wave  is  formed  by 
the  corner  deflection  which  interacts  with  the  turbu¬ 
lent  boundary  layer.  At  sufficiently  large  corner  angle, 
the  boundary  layer  separates  at  the  corner  and  shock 
bifurcates  into  separation  and  reattachment  compres¬ 
sion  systems.  The  test  conditions  for  the  selected  cases 
are  shown  in  Table  13.  Experimental  data  of  Zhel¬ 
tovodov  et  a/®®“®^  includes  surface  mean  pressure,  skin 
friction,  adiabatic  wall  temperature  and  heat  transfer; 
and  boundary  layer  profiles  of  mean  velocity  and  tem¬ 
perature  and  turbulence  statistics. 


Table  13  Supersonic  Compression  Corner 

Quantity  Value 

Moo  2.9 

Res  7.5  X  10^  to  1.5  x  10® 
a  8°  and  25° 

Dataset  No.  2:  2-D  Expansion-Compression 
Corner 

Prepared  by 
A.  Zheltovodov 

Description  of  Experiments 

A  two  dimensional  equilibrium  supersonic  turbulent 
boundary  layer  develops  on  a  flat  plate.  It  expands 
through  an  angle  a  at  a  corner  and  subsequently  is 
compressed  by  an  equal  angle  as  illustrated  in  Eig.  19. 
The  expansion  reduces  the  turbulence  kinetic  energy 
and  increases  the  mean  kinetic  energy  in  the  boundary 
layer.  The  turbulence  structure  within  the  boundary 
layer  is  no  longer  in  equilibrium  upon  reaching  the 
compression  corner.  Depending  on  the  compression 
corner  angle,  Mach  number  and  Reynolds  number,  the 
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boundary  layer  may  separate  at  the  compression.  Ex¬ 
perimental  data  of  Zheltovodov  et  is  available. 
The  experimental  condition  for  the  selected  case  are 
shown  in  Table  14.  Additional  experimental  data  is 
also  available. 


Table  14  Expansion  Compression  Corner 
Quantity  Value 

2.9 

Res  4.1  X  10^  to  1.95  x  10® 
a  25° 


Dataset  No.  3:  2-D  Shock  Impingement 

Prepared  by 

J.  Debieve,  M.  Elena,  P.  Dupont 

Description  of  Experiments 

A  two-dimensional  supersonic  equilibrium  turbulent 
boundary  layer  develops  on  a  flat  plate.  It  is  subjected 
to  a  shock  wave  produced  by  a  shock  generator  (wedge) 
of  angle  a  placed  in  the  external  flow,  as  indicated 
in  Eig.  20.  Eor  sufficiently  large  angle  of  yaw  of  the 
generator,  separation  occurs.  Although  the  incident 
shock  is  steady,  the  reflected  shock  shows  significant 
oscillations.  The  boundary  layer  is  fully  turbulent,  but 
the  Reynolds  number  is  rather  low,  so  that  the  data  are 
appropriate  for  comparisons  with  LES  computations. 


The  fiow  conditions  are  listed  in  Table  15.  The 
incoming  boundary  layer  is  fully  turbulent  and  devel¬ 
ops  on  a  fiat  plate  with  nearly  adiabatic  constant  wall 
temperature.  Experimental  data  are  available  for  gen¬ 
erator  angles  a  =  4°  and  8°.  Experimental  data®®^®® 
include  surface  mean  pressure  and  mean  temperature, 
boundary  layer  profiles  of  mean  velocity  and  temper¬ 
ature,  three  Reynolds  stresses  (turbulent  friction  and 
two  normal  stresses) ,  and  turbulence  statistics  on  tem¬ 
perature  fiuctuations.  Conditional  statistics  have  been 
performed  to  educe  the  large  scale  eddies  developed 


downstream  of  separation.  Statistics  of  the  Strouhal 
number  related  to  these  large  scale  structures  (pdf) 
have  been  obtained. 

Table  15  Shock  impingement 

Quantity  Value 

Moo  2.3 

Res  5.2  X  10^ 

a  4°  and  8° 


Dataset  No.  4:  3-D  Single  Fin 

Prepared  by 
D.  Knight 

Description  of  Experiments 

A  two  dimensional  supersonic  equilibrium  bound¬ 
ary  layer  develops  on  a  fiat  plate.  The  deflection  of 
the  fiow  by  a  wedge  of  angle  a  (Eig.  21)  generates 
an  oblique  shock  which  interacts  with  the  turbulent 
boundary  layer.  An  extensive  set  of  experiments  have 
been  performed  for  this  configuration.  Three  exper¬ 
iments  were  selected  for  CED  validation  (Table  16). 
The  experimental  data  for  Case  Nos.  1  and  3  are  from 
Zheltovodov  et  al,^^’^^  and  the  data  for  Case  No.  2  is 
from  Kim  et  al.^^  Experimental  data  include  surface 
pressure,  skin  friction  and  surface  streamlines. 


Fig.  21  3-D  single  fin 


Table  16  3-D  Single  Fin 


Case 

Moo 

a 

Res 

1 

3.0 

15° 

1.9  X  10® 

2 

4.0 

20° 

2.1  X  10® 

3 

4.0 

30.6° 

1.6  X  10® 

Dataset  No.  5:  3-D  Double  Fin  Interaction  No.  1 

Prepared  by 
D.  Knight 

Description  of  Experiments 

A  two-dimensional  supersonic  equilibrium  turbulent 
boundary  layer  flows  into  a  rectangular  channel  formed 
by  two  fins  of  angles  cei  and  a2  attached  normal  to 
the  fiat  plate  as  shown  in  Eig.  22.  Shock  waves  are 


formed  by  the  fin  deflections  which  interact  with  the 
turbulent  boundary  layer  on  the  flat  plate.  For  suf¬ 
ficiently  strong  shocks,  the  boundary  layer  separates 
into  counter-rotating  vortices  which  intersect  to  form 
a  vortex  pair.  A  complex  pattern  of  shock  waves,  ex¬ 
pansions  and  slip  lines  form  due  to  the  intersection  of 
the  incident  shocks.  A  low  total  pressure  region  coin¬ 
cides  with  the  vortex  pair.  Secondary  separation  may 
occur  downstream  of  the  intersection  of  the  primary 
separation  lines. 


Fig.  22  3-D  double  fin 


The  experiments  were  performed  at  the  Penn  State 
Gas  Dynamics  Laboratory  by  Settles  et  The 

incoming  boundary  layer  is  an  equilibrium  turbulent 
boundary  layer  with  an  adiabatic  or  isothermal  wall. 
The  flow  conditions  are  indicated  in  Table  17.  Experi¬ 
mental  data  includes  surface  pressure,  surface  flow  vi¬ 
sualization,  Laser  Interferometer  Skin  Friction  (LISE) 
measurements  of  surface  shear  stress,  flowfleld  pitot 
pressure  surveys,  and  Planar  Laser  Scattering  (PLS) 
images  of  the  flowfleld. 

Table  17  3-D  Double  Fin 


Case 

-^oo 

Qi 

02 

Res 

1 

3.9 

9° 

9° 

2.6  X  10 

2 

3.9 

11° 

11° 

2.6  X  10 

3 

3.9 

13° 

13° 

2.6  X  10 

4 

3.9 

15° 

15° 

2.6  X  10 

Dataset  No.  6:  3-D  Double  Fin  Interaction  No.  2 

Prepared  by 

A.  Zheltovodov,  A.  Maksimov. 

Description  of  Experiments 

The  experiments  were  performed  at  the  Institute  of 
Theoretical  and  Applied  Mechanics  by  Zheltovodov  et 
The  incoming  boundary  layer  is  an  equilib¬ 
rium  turbulent  boundary  layer  with  an  adiabatic  or 
isothermal  wall.  The  flow  conditions  of  the  selected 
experiments  are  given  in  Table  18.  Experimental  data 
is  also  available  for  other  configurations.  Experimental 
data  includes  surface  pressure,  adiabatic  wall  temper¬ 
ature,  heat  transfer  and  surface  flow  visualization. 


Table  18  3-D  Double  Fin 


Case 

-A^oo 

ai 

02 

Res 

1 

4.0 

7° 

11° 

3.0  X  10' 

2 

4.0 

15° 

15° 

3.0  X  10' 

Dataset  No.  7:  3-D  Double  Fin  Interaction  No.  3 

Prepared  by 

A.  Zheltovodov,  E.Schiilein. 

Description  of  Experiments 

The  experiments^^^^®  were  performed  at  the  DLR 
Institute  of  Fluid  Mechanics,  Gottingen,  Germany. 
The  configurations  are  shown  in  Figs.  23  and  24  where 
the  dimensions  are  in  mm.  The  incoming  boundary 
layer  is  an  equilibrium  turbulent  boundary  layer  with 
an  isothermal  wall  (T^j  =  290  —  300K)  at  Mach  5, 
Res  =  (1.9  —  2)  X  10®,  T^fTt  =  0.691.  Experimental 
data  is  available  for  fins  angles  fi  =  18°  and  23°  (Ta¬ 
ble  19).  Experimental  data  includes  surface  pressure, 
surface  flow  visualization  and  optical  crossing  shock 
wave  structure  visualization. 


Fig.  24  Long  fins  (/3  =  23°) 


Table  19  Double  Fin 


/? 

Type 

L 

A 

B 

1 

8° 

SF 

229.5 

142.0 

100 

309 

12° 

SF 

249.5 

161.5 

100 

280 

18° 

SF 

269.0 

196.0 

100 

260 

23° 

LF 

286.0 

307.0 

100 

dimensions  in  mm 


Dataset  No.  8:  Large  Cone-Flare 

Prepared  by 

M.  S.  Holden 

Introduction 

An  experimental  study^®’ in  which  surface  and 
flowfield  measurements  were  made  has  been  conducted 
to  examine  the  structure  of  turbulent  flow  separation 
over  large  cone-flare  configurations.  This  study  was 
conducted  in  Calspan’s  96  inch  Shock  Tunnel  at  Mach 
numbers  of  11  and  13  with  Reynolds  numbers  up  to 
100x10®.  The  measurements  made  at  the  large  unit 
Reynolds  number  and  models  employed  in  this  study 
demonstrated  that  the  attached  and  separated  flows 
were  fully  turbulent  and  a  thick  boundary  layer  was 
developed  at  the  cone/flare  junction,  which  enabled 
highly  resolved  measurements  in  the  interaction  re¬ 
gion.  Surface  heat  transfer  and  pressure  measurements 
were  made  in  attached  and  separated  flows  at  the 
cone-flare  junction  for  30°  and  36°  flare  angles  respec¬ 
tively.  Flowfield  surveys  were  made  in  the  separation 
region  with  pitot  pressure  and  total  temperature  rakes. 
Holographic  interferometry  and  Schlieren  photography 
were  used  to  obtain  details  of  the  flowfield  structure. 
This  study  suggests  that,  in  hypersonic  flow,  the  sep¬ 
aration  region  extends  only  a  very  small  fraction  of 
the  boundary  layer  thickness  and  is  a  highly  unsteady 
process.  Only  by  employing  instrumentation  with  fre¬ 
quency  response  fast  enough  to  follow  the  unsteady 
movement  of  the  separation  shock  is  it  possible  to  de¬ 
termine  the  fundamental  structure  of  flow  separation 
in  turbulent  hypersonic  flows. 

Description  of  Experiments 

For  this  study,  a  large  6°  cone  with  flares  of  ei¬ 
ther  30°  or  36°  attached  to  its  base  were  employed 
in  the  experiment.  The  cone  angle  and  model  length 
were  selected,  on  the  basis  of  calculations,  to  achieve 
the  maximum  length  over  which  uniform  constant- 
pressure  flow  could  be  established  within  the  tunnel. 
The  model  configuration  is  shown  in  Fig.  25. 


The  flowfield  is  a  perfect  gas,  axisymmetric  flowfield, 
with  a  well  developed  turbulent  boundary  layer.  Shock 
wave-boundary  layer  interaction  and  flow  separation 
constitute  the  principal  flowfield  physics. 

This  dataset  contains  measurements  at  the  largest 
Mach  number  and  Reynolds  number  currently  avail¬ 
able.  The  experimental  data  includes  surface  temper¬ 
atures,  heat  transfer  and  pressure  distributions.  Addi¬ 
tionally,  the  use  of  a  single-pass  Schlieren  system  with 
a  focal  length  of  10  feet  for  flow  visualization  was  used. 
The  freestream  conditions  are  shown  in  Table  20. 


Table  20  Large  Cone-Flare 


Runs 

4,8 

6,7 

Mi 

3.345 

3.633 

Po  (psia) 

7.216  X  10^ 

1.760  X  10^ 

Hq  (ft^/sec^) 

1.825  X  10^ 

2.147  X  10^ 

To  (°R) 

2.717  X  10^ 

3.104  X  10® 

M 

10.96 

13.01 

U  (ft /sec) 

5.922  X  10^ 

6.458  X  10® 

T  (°R) 

121.4 

102.6 

Poo  (psia) 

9.172  X  10-2 

7.345  X  10-2 

q-oo  (psia) 

7.721 

8.712 

Poo  (slug/ft®) 

6.340  X  10“® 

6.038  X  10-® 

Poo  (slug/ft-sec) 

1.021  X  10“'^ 

8.634  X  10-® 

Re  (ft-i) 

3.680  X  10® 

4.544  X  10® 

Po  (pitot)  (psia) 

1.431  X  10^ 

1.619  X  10^ 

6.  BASE  FLOWS  WITH  AND  WITHOUT 
PLUME  INTERACTIONS 
TEAM  LEADERS:  P.  BARKER 
P.  REIJASSE 

Prepared  by 

P.  G.  Bakker  and  P.  Reijasse 
Introduction 

For  the  development  of  the  next  generation  of 
reusable  launchers  and  reentry  vehicles,  one  of  the  crit¬ 
ical  areas  is  the  proper  modelling  of  the  base  flow.  The 
low  pressures  that  act  on  the  base  region  of  bodies  in 
supersonic  and  hypersonic  flight  can  cause  significant 
amounts  of  drag.  Due  to  separation  and  reattachment 
in  the  baseflow  region  also  the  heat  loading  at  the  base 
may  be  considerable.  Futhermore  non-axisymmetric 
effects  and  unsteadiness  can  cause  severe  sideloads  on 
a  launch  vehicle  in  the  ascent  phase.  Important  pa¬ 
rameters  influencing  the  base  flow  physics  are,  e.g., 
boattailing  of  the  afterbody,  base  bleed,  presence  of  ex¬ 
ternal  flow  and  conditions  of  the  exhausting  jet-plume 
(underexpanded  or  overexpanded).  Three  different 
experimental  data  sets  meeting  the  WGIO  objectives 
were  identified  covering  most  of  these  effects  and  (par¬ 
tially)  their  interactions.  Dataset  No.  1  (TU-Delft) 
is  a  cylindrical  afterbody  with  jet  in  external  flow. 
Dataset  No.  2  (Univ.  Illinois)  is  a  cylindrical  and  boat- 
tailed  afterbody  in  external  flow  with  no  jet.  Dataset 
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No. 3  (ONERA)  is  a  boattailed  afterbody  with  jet  in 
external  flow.  Herrin  and  Dutton  (Dataset  No.  2) 
studied  the  influence  of  boattailing  in  absence  of  a 
jet-plume.  Their  experimental  results  show  that  boat¬ 
tailing  (without  plume)  raises  the  base  pressure  level 
and  lowers  the  turbulence  levels  in  the  shear  layer,  so 
a  signiflcant  reduction  in  net  afterbody  drag  results. 
Dataset  Nos.  1  and  3  both  deal  with  external  flow 
and  an  underexpanded  jet-plume.  Comparing  them 
we  can  study  the  influence  of  boattailing  in  presence 
of  a  jet-plume  (underexpanded). 

Dataset  No.  1:  Base  Flow-Plume  Interaction  in 
Supersonic  External  Flow 

Prepared  by 

P.  Bakker  and  W.  Bannink 

Description  of  Experiments 

The  test  setup  enables  the  investigation  of  the  flow- 
fleld  along  an  axisymmetric  body  with  a  single  operat¬ 
ing  exhaust  nozzle.  The  supersonic  jet  emanating  from 
the  centrally  protruding  exhaust  nozzle  in  the  base  in¬ 
teracts  with  the  external  main  flow.  Essentially  two 
independant  supersonic  streams  -  the  external  flow 
and  the  central  jet  -  interfere  at  the  base  region.  In  the 
interaction  zone  a  turbulent  mixing  layer,  a  recirculat¬ 
ing  region  and  a  shock  system  (plume  shock,  barrel 
shock,  Mach  disc)  is  formed.  The  model  has  a  conical 
forebody  (semi-apex  angle  of  11°)  with  a  fair  amount 
of  bluntness  (nose  radius  of  7.5  mm)  and  a  cylindrical 
afterbody  (diameter  of  50  mm,  length  of  90  mm) .  The 
total  length  of  the  model  is  186.81  mm.  The  model  is 
strut-mounted  at  the  lower  side  and  has  a  free  base. 
Erom  the  centre  of  the  base  a  nozzle  protrudes;  its  out¬ 
side  shape  is  a  circular  cylinder.  The  nozzle  itself  is 
conical  (total  divergence  of  15°)  having  Mach  4  at  the 
exit.  The  model  is  positioned  in  a  uniform  supersonic 
free  stream  at  zero  incidence.  Additional  details  are 
presented  in  Bannink  et 

The  experimental  conditions  are  shown  in  Table  21. 
The  boundary  layer  along  the  body  is  tripped  at  38 
mm  from  the  nose.  Experiments  are  performed  at 
jet  pressures  (total  pressures)  ranging  from  1.5  MPa 
to  10  MPa.  Due  to  the  unexpected  low  base  pres¬ 
sures  it  turned  out  that  this  pressure  range  covers  only 
underexpanded  jets.  Apart  from  these  experiments, 
experiments  without  jet  are  performed  as  well.  Exper¬ 
imental  data  include  mean  static  surface  pressures  on 
afterbody  and  on  the  base,  Schlieren  and  shadowgraph 
pictures  of  the  base  flow,  and  mean  Pitot  pressures  in 
the  base  area  to  capture  shock-system  and  shear  lay¬ 
ers.  New  experiments  both  for  Mqo  =  2  and  Mqo  =  3 
are  underway  at  jet  pressures  ranging  from  0.3  to  1.2 
MPa  in  order  to  capture  also  overexpanded  exhaust. 
Expected  experimental  data  will  be  similar  to  those 
available  for  the  underexpanded  case. 

Some  computational  results  of  two  codes  based  on 


Table  21  Base  Flow  Plume  Interaction 


^oo 

Rcl 

Pt^ 

Tt^ 

2 

5.16  X  10® 

0.206  MPa 

285°K 

3 

8.70  X  10® 

0.575  MPa 

285  K 

RANS 

equations  are 

available. 

Houtman  (TU- 

Delft)  has  performed  numerical  simulations  using  a 
code  based  on  a  cell  centered  finite  volume  discretiza¬ 
tion  of  the  Euler-  and  the  Navier-Stokes  equations. 
The  code  is  equipped  with  the  one-equation  Spalart- 
Almaras  turbulence  model.  Structured  multiblock 
grids  are  used.  Van  der  Weide  (von  Karman  Insti¬ 
tute)  has  used  a  code  based  on  a  multidimensinal 
upwinding  method  on  unstructured  grids.  Various  tur¬ 
bulence  models  {e.g.,  Spalart-Almaras,  BSL  and  SST) 
are  used.  The  computational  results  only  show  good 
agreement  with  experiments  in  the  not  viscous  domi¬ 
nated  parts  of  the  flowfleld.  The  prediction  of  the  base 
flow  however,  appears  a  tough  problem  which  requires 
extremely  fine  grids,  especially  for  the  two-equation 
models.  The  calculated  base  pressures  are  still  signifi¬ 
cantly  lower  than  those  found  in  the  experiments. 

Dataset  No.  2:  Supersonic  Axisymmetric  Base 
Flows  Including  Boattail  Effects 

Prepared  by 

P.  Bakker  and  W.  Bannink 

Description  of  Experiments 

The  configuration  enables  detailed  non-intrusive 
measurements  on  the  entire  near-wake  flowfleld  struc¬ 
ture  behind  a  cylindrical  and  a  boattailed  afterbody 
immersed  in  a  supersonic  flow.  The  near-wake  is  char¬ 
acterized  by  the  occurrence  of  separating  shear  layers 
reattaching  downstream,  and  enclosing  a  recirculation 
region.  The  base  model  is  mounted  on  a  central  sting 
which  supports  the  model  from  upstream.  The  cylin¬ 
drical  part  of  the  afterbody  used  in  the  experiments  is 
63.5  mm  diameter. 

Table  22  Base  Flow 

Moo  Re/m  Pt^  Tt^ 

2.46  5.2  X  10^  0.515  MPa  294  K 

Experiments  were  performed  at  the  University  of 
Illinois  Urbana-Champaign  by  Dutton  et  af.*®  **  The 
flow  conditions  are  shown  in  Table  22.  Oilstreak  vi¬ 
sualization  and  micrometer  measurements  were  per¬ 
formed  to  attain  an  axisymmetric  flow.  Experimental 
data  include  Schlieren  and  shadowgraph  pictures  of 
the  near-wake  flow,  mean  static  pressures  at  several  lo¬ 
cations  on  the  base  and  on  the  afterbody  surfaces,  and 
two-component  LDV  measurements  in  the  boundary 
layer  approaching  the  base  corner  and  in  the  near-wake 
region.  Erom  the  detailed  quantitative  data  obtained 
in  the  near-wake,  specifically  afterbody  and  base  pres¬ 
sure  distributions,  mean  velocities,  turbulence  intensi¬ 
ties  and  Reynolds  shear  stresses  have  been  obtained. 
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Dataset  No.  3:  Plume-Induced  Flow  Separation 
on  a  Jet-On  Axisymmetric  Boattailed  Afterbody 

Prepared  by 

P.  Reijasse,  B.  Corbel 

Introduction 

These  experiments*®  have  been  performed  under 
contract  of  the  French  Ministry  of  Defense  to  consti¬ 
tute  a  Navier-Stokes  validation  test  case  regarding  the 
pluming  phenomena  of  an  underexpanded  propulsive 
jet  which  induces  the  external  flow  separation  on  a  ax¬ 
isymmetric  afterbody  boattail.  The  tests  have  been 
carried  out  in  the  S8Ch  continuous  atmospheric  sub¬ 
sonic/supersonic  research  windtunnel  in  the  Funda¬ 
mental  and  Experimental  Aerodynamics  Department 
of  ONERA,  headed  by  Prof.  J.  Delery. 

Description  of  Experiments 

The  S8Ch  wind  tunnel  was  equipped  with  two  pla¬ 
nar  converging-diverging  half- nozzles.  The  test  sec¬ 
tion  is  0.12  m  X  0.12  m.  The  axisymmetric  model  is 
supported  by  an  upstream  central  sting  fixed  in  the 
settling  chamber.  A  sketch  of  the  model  is  shown  in 
Eig.  26.  The  external  diameter  D  of  the  model  is  30 
mm.  The  aft  part  of  the  model  is  a  9B05  boattail 
having  a  length  of  31.5  mm.  The  model  is  equipped 
of  a  lOBO  conical  nozzle  having  an  exit  diameter  of 
14.9  mm.  The  origin  of  the  coordinate  system  is  lo¬ 
cated  in  the  exit  plane  of  the  model  nozzle  on  the 
symmetry  axis.  The  mesh  coordinates  are  nondimen- 
sionalized  by  the  model  diameter  D. 


Fig.  26  Jet-on  boattailed  afterbody 

The  afterbody  region  of  the  model  is  subject  to  com¬ 
plex  interactions  due  to  flow  confluence  between  the 
propulsive  jet  and  the  outer  flow.  After  the  upstream 
external  flow  crosses  the  expansion  waves  centered  on 
the  boattail  origin,  it  meets  an  interaction  shock  lo¬ 
cated  at  the  middle  of  the  boattail.  This  shock  wave 
/  boundary  layer  interaction,  resulting  in  the  separa¬ 
tion  of  the  external  boundary  layer,  is  induced  by  the 
pluming  of  the  jet,  with  the  underexpanded  jet  behav¬ 


ing  as  a  “fluid  ramp”  facing  to  the  incoming  external 
flow.  The  interaction  shock  coalesces  with  the  conflu¬ 
ence  shock  crossing  the  external  flow.  Downstream  of 
the  nozzle  exit,  the  jet  is  accelerated  when  it  crosses 
the  expansion  waves  centered  on  the  nozzle  lip.  The 
thin  jet  shear  layer  rapidly  thickens  when  it  reaches 
the  confluence  region.  Then  the  confluence  process 
starts  the  development  of  a  wake  which  ensures  the 
mutual  adaptation  of  the  flows.  The  jet  internal  shock 
-  or  barrel  shock  -  which  results  from  the  Mach  lines 
focusing  process,  reflects  on  the  symmetry  axis  in  a 
singular  way  with  the  formation  of  a  Mach  disk  struc¬ 
ture. 

The  experimental  database  includes  wall  pressure 
on  the  afterbody  boattail  and  on  the  base,  spark  and 
continuous  Schlieren  photographs,  and  2000  measure¬ 
ment  points  by  the  two-component  laser  velocimeter 
technique  giving  the  mean  streamwise  and  radial  com¬ 
ponents  of  the  velocity,  the  streamwise  and  radial  tur¬ 
bulence  intensities,  the  turbulent  shear  stress  and  an 
approximation  of  the  turbulence  kinetic  energy. 

Ereestream  conditions  are  given  in  Table  23.  The 
unit  Reynolds  number  is  equal  to  12.24  x  10®  m^^. 


Table  23  Jet-on  Boattailed  Afterbody 


Flow  parameter 

External  Flow 

Jet  Flow 

Specific  heat  ratio 

1.4 

1.4 

Pt  (bar) 

0.975 

7.75 

Tt  (K) 

298 

298 

M 

1.94 

1.75 
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